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Dr. Dennis Bogan 

NASA, Office of Space Science 
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Two Independence Square 

Washington, D.C. 20546-0001 


RE: Renewal of NASA NAGW-4388, January, 1997-December, 1997 


Dear Jay and Dennis: 

Much of my time during the first five months of 1996 were perforce devoted to 
analyzing data acquired by the Galileo Probe Mass Spectrometer (GPMS) and writing the 
GPMS team's Science paper. Because of my 18 years of experience analyzing mass 
spectrometer data from the Pioneer Venus Probe Mission (LNMS), and my emeritus status 
I found myself uniquely prepared to make rapid progress in analyzing and interpreting the 
GPMS results and with time to do so. I enclose a copy of the Science paper that resulted. 
Since that paper was written, post flight calibration is proceeding methodically and we have 
had a careful look at the data obtained between 14.5 bar and 21 bar which was barely 
discussed in the Science paper. (See the COSPAR paper, also enclosed) The most 
interesting result here is that in this altitude range the water vapor mixing ratio gets up to 
1.6 times solar, compared to the 0.2 times solar below 11 bar that was quoted in our 
preliminary report. Like H 2 S and HC1, H 2 O increases in mixing ratio with decreasing 
altitude (or perhaps horizontal displacement). We are inclined to attribute this effect to the 
probe’s finding itself in the descending volatile poor arm of a circulation system. We are 
awaiting calibration of the instrument for NH 3 for further enlightenment. 

In the meanwhile, I have been able to get started on a couple of Venus hydrogen 
projects. On of these should be finished and on the way to publication by the time this 
letter would have been written in the past. I have been trying to understand how to 
reconcile the standard treatment of the evolution of the H 2 O and HDO reservoirs on the 
planet over 4.5 Gyr in the presence of H and D escape and injection by comets. The simple 
relations 
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govern this evolution. 

Where the fractionation factor f is given by 

f MSI 

4>l/[H] ' 

and, if we take 

4>1=-K[H] , 

with K constant (big assumption^and let 

Rs = (D/H), 

in the cometary source, the solutions to (1) and (2) are 

[Hl/[H 0 ] = e-^ + £l(M^(,. e -K.) 


(3) 

(4) 

(5) 

(6) 
(7) 


[D] (l-e-^j/f + (R s /R 0 XP 1 /»im/H o ])~ 1 e- fKt 

[H] 1 ~ e_Kt + (Pi/ 1 <(>! [H]/[H 0 ]) -1 e _Kt 

These relations are discussed in detail in a chapter for Venus II I have written (copy 
enclosed). 

Early this year my GSFC PV colleagues and I published a paper in JGR, Planets 
(Hartle et al, enclosed) showing how analysis of PVO data over several solar cycles 
allowed us to determine H and D escape fluxes <J>i and <j >2 averaged over the planet and the 
solar cycle: 

01 = 7xl0 6 cm -2 s -1 , 
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f = 0.44. 

This is close to the Charge Separation Electric Field H + flux obtained by Hartle and 
Grebowski and generalized by Hartle and myself for solar cycle effects 

<J >1 (E) = 9xl0 6 cm -2 s -1 

f = 0.17 

It is not a good representation of my "best guess" for the combined electric field, charge 
exchange flux (E + CE) (See Hartle et al, Donahue and Hartle, Donahue et al, all enclosed). 

<t>i(E + CE) = 16xl0 6 cm _2 s _1 , 

f = 0.10. 

We can put these values for <bj and f in (6) and (8), along with the measured value of 

R = D/H = 150 ±30 
and 

P! = 7xl0 5 - 4.2xl0 6 cm -2 s -1 


(see Donahue et al, Venus II), 
where 


1<R S /R 0 <2 


The funny thing is that the problem is over determined if we have [H], <pj , R, f, 
Pi, Rs/Ro given by experiment. This is because the measurement of R at 4.5 Gyr and 

K = {0>i/[H]} 45 


means, for any assumed Rs/Ro 


ZlM 

fl [Ho] 


'4.5 


is determined. Then in (6) 
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and [H 0 /H] which gives the ratio of the original water inventory to the present on is set. 
There is no freedom in setting Pj/^x ! 

I am now in the process of investigating the implications. 

I am also preparing to have a go at calculating the charge exchange contribution to 
hydrogen loss rates, using realistic models for exospheric H, H + , D, D + , and ion 
temperature from PV data. This has never been done properly. Hodges and I have 
proposed doing this a couple of times, but Hodges has not been funded for his part. I hope 
to get his "unfunded" help in closing this embarrassing gap in our knowledge. 

I have presented papers on Cytherean and Martian Water at the European 
Geological Society Meeting at the Hague and at the Pale Blue Dot workshop recently at 
NASA ARC. 

During the final year of this grant I propose: 

a) to finish the calculation of charge exchange loss rates of H and D from Venus; 

b) to continue to analyze Galileo Probe Neutral Mass Spectrometer data after 
MODA support terminates on April 30, 1997. 

c) To begin work on the sulfur chemistry and redox state of the atmosphere of 
Venus. 
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The Galileo Probe Maas Spectrometer: 
Composition of Jupiter’s Atmosphere 

Hasso B. Niamann, Suahil K. Atraya, Qaorje R. Carignan, 
Thomas M. Donahue, John A. Habsrman, Dan N. Harpold, 
Richard E. Hartle, Donald M. Human, Wayne T. Kaaprzak, 
Paul R. Mahaffy, Tobias C. Owen. Nelson W. Spencer, 
Stanley H. Way 

Tha composition of tha Jovian atmoaphara from in altitude markad by in atmosphartc 
prtaaura of 0.3 to 21 bars waa dtiarmlnad by a quadrupola maaa apactromatar on tha 
Gawao prooa. Tna mixing r«no of Ht (naiium) to H a (hydrogan), o.i 50, la eioaa to tna aoiar 
ratio. Tha aounoanoaa of mathana, watar, argon, naon, and dlhydrogan aulfida war* 
maaaurad; Krypton and xanon wart datactad. Aa maaaurad In tha Jovian atmoephare, tha 
amount of carton a 8.9 timaa tha aoiar abundanea relative to H a , tha amount of aulfur 
la groa t arthan tha aoiar atoundanoo, and tha amount of oxygon la mueh loot than tha aoiar 
abundanea. Tha naon abundanea eomparad with that of hydrogan la about an erdar of 
magnltuda lau than tha aoiar abundanea. Isotopic ratios of carbon and tha nobla gaaaa 
art conaiatant with aoiar valuaa. Tha maaaurad ratio of dautartum to hydrogan (0/H) of 
(S ± 2) x 10~* Indieataa that this ratio la graatar In sotar-syatam hydrogan than In local 
Intaratalar hydrogan, and tha *HaTHa ratio of (1 .1 ± 0.2) x 10~* provldaa a now vaiua 
for protoaoiar (aoiar naouia] haiium laotopaa. Togatnar, tna D/M and “NaTHa ratios ara 
conalatart with oonvaralon In tha sun of protoaoiar dautartum to praaant-day *Ha. 


Determination of the compmlrinn of Junl- 
ter’s atmosphere ihould coru train the rel»- 
dve Importance of direct contribution* to 
tha atmotphata from tha tolar nebula itself, 
on tha one hand, and from largt Icy or 
rocky object' present in the early outtr soUr 
system, on the other. The contribution of 
(hate "planeteiimaU" could have been In 
tha form of an early stmo*ph*r» around 
Jupiter’s primitive lolid com or in tha form 
of vnlatllaa they carried to the planet after it 
began to acquire la atmotphata of H and 
He from the tolar nebula (1,2). Tha dagraa 
of resemblance between the lun'i atmo- 
tphate and that of Jupiter in tha abundanc- 
e* relative to H of alamano such ai C, N, 
and O ihould be dectilve in providing the 
inquired corocralnti. 

Before the direet. In litu m#**ursmant» 
reported here, remote tpeetroscoplc renting 
from Barth and from spacecraft had indicate 
td that the relative abundance! of C (in 
CH 4 ) and N (In NH.) wan gnatar than 
the aoiar abundances (3, 4). Phorphlna and 
water vapor had been detected, but their 
abundance! below tha cloud* were uncer- 
tain, and H,S and die noble gaeci other 
than He had not been detected at ell, De- 
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Mongan, 2«M iieywwu strew. Am Aruvr. mi aeioa, 
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determining the iiotopic ratio) of their con- 
itituent eltmeno (5). Signal* ware recorded 
from over 6000 individual value* of mm*. 
to-charg« ratio [nJt). The data returned by 
the Instrument have led to tha ditcovety of 
•ix new atuuMplieiic unutitwnn and the 
mceturemcnt of numcroui abundance* and 
Itotopic retie*. Tha gas sampling tyitem 
admirtad rhe Jovian armnaphrre through 
low-conductance leak) to the ion source of 
a quadrupole mail analyser. Tha Ion lourea 
wai pumped by a gttter and the anilyw 
volume pumped to a much lower proeure 
by a getter backed by a iputttr-ion pump. 
One tnlec (direct teak 1; DL1) waa open 
from 0,22 to 2.76 ban; the odist Inlet 
(DL2) functioned from 8.21 ben to the end 
at about 2 1 ben ( 6 ). 

Figure 1 shows a temple mass ipectnun 
from tha 3-bat region at 228 K, before any 
•ignificant water cloud wa* expected, and 
another one from the 11-bar region, where 
the 390 K temperature axurti that the 
probe waa far below any watar condamacion 

level. Cumpariaun uf die two spectra shows 
a large increaac in the iccond direct leek of 
HjO end NH). In fact, the detign of the 
leak 1 ump ling system, with a relatively 
flnltive data to replace theu tantaliiing long vacuum path to tha Ion tource, pre- 
duet to the taurcei of volatile compcundi -• elude* a aeniltive detection of the surface 
In Juplter^tmoiphere were not ayatfeH^tvwctivs ipeciei H 2 0 and NH, In the atreo- 
before tbt^deicent cf firtjptahrr -fjjt . .«pherlc region wrapled by this leak. Raw 

An accurate direct meaiurement of the count* ihown in there ipectra must be cor- 
abundance* of He and the other noble garei rected for the counting tyitcm dead time. 

It important for undentanding how Jupiter The efficiencies of die ion tource and the 
hat proctucd it) constituents. Abundance* vpumping tyitcm vary from ipeclet to epe- 
law than tha n of die tolar vaiuci would/ cle», and the icnaitivity factor*, which ate 
tignal the formation, in die metallic inwri- different for the rwn direct leak), mutt be 
or of the planar, of a reparate He phase, 
tome of which hot precipitated deeper into 
the interior. Remote tensing had indicated 
a modeit depletion of He, but this Indirect 
result is much lei* satisfactory than the 
direct tensing by two of tha Galileo Instru- 
ment) repotted here. 

The Galileo probe mat* ipcctromster 
(OPMS) wu deiigned to measure the mix- 
ing ratio) of major and minor tpecies while 

TaDM 1, MNrg ratio* to H t and isotopic ratio* In to* jovien Mr: wspi wre. bier* apece* mdlcet* *w. 

Ratio to 


individually determined. In most cases, we 
used preflight calibrations, but for die moat 
important ipacltt we were able to do some 
recalibrations on the spare Instrument. 
With the direct leaks, the background level 
was typically between 0 and 3 count*. Two 
background measurement) demonstrated 
that the instrument arrived at Jupiter in a 
very clean ttatt end quickly Inst it* memory 
of almost all gate) contributed by DL1. 
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Fla. 1.5tmpltnraa$ptctraTr>«ipictmln(A|Wiftodt»i.naotfvouo^C)Ul Cftw*anpfawuntof2.72 ■rda.OSbar. wttn wnofrrrt t*mp*r»iur« ti 216 K. Tnm 
•ptdra In (I) ww ootaintd through OH batwaan pmturtt of i 0,6** 1 1 .3 bw, with tmbiftrt tompmium fit 360 K. No oorreotlorw tmabiin ratter *id 
«m* or for dW teftn t tfWonoiot «t dftef*nt mooooo. 
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TabU 1 ihowt th* mixing ratio* (num- 
ber densities relative to H,) of ipedei (hat 
were detected and for which tellable mea- 
surement. could be made. Error csdmatcs 
•re roughly ltr and reprewnt the data scat- 
ter about the average value. Some of the 
obaerved mixing ratioi can be compared 
with three for the iwler .y.tem a* compiled 
by Anders and Gtcvcmc (7), commonly 
referred to a* the tolar abundant.. Abo 
thown ate multi nf earlier analyses baaed 
on Voyager or Earth-based ipectra. 

The jovian H«/H, ratio, 0.156 - 0.006, 
it die avenge of 14 Individual measure- 
menu from DLl (0.157 ± 0.006) and 4 
from DL2 (0.151 t 0.004). All die other 
noble gaaei were detecred. and they .re 
elementally fractionated. The Ne ratio to 
A. 1» about an order of magnitude leu then the 
tolot ratio, whotaaa for Ar the ratio i* .light- 
ly greater. Kr end Xe warn dearly detected: 
upper limita on the mixing ratioi ere teepee* 
lively 5 end 50 time, the tolar value*. For 
all four of die heavier noble taees, die 
iaotopic ratioi are dote to the aolar valuer. 

Tha ’HefHc and the D/H ratioi were 
determined from an analyala of the 3 «nu 
(atuinic mau unit) data. This analyri. ir 
complicated by tha need to correct for Hj + 
produced in the ion mure, end to disfin- 
between the contribution* of HD end 
He. H,+ la generated in ana of two wayi; 
by dissociative ionisation of CH 4 or by • 
two-body ion-molecule reaction between 
Hj + end H r The efficiency of produc- 
tion could be deterrained during the second 
backgtuund study whan little Hi and He, 
but a contld.tabi. amount of CH,, «xiired 
in tha ion Kairca and CH 4 wat the only 
important lource of 3 -emu ion*. During pre- 
K flight calibration, when the Hcpreoure in 
& the CPMS wa* varied over a'lfcge range. 
k the HD/H, ratio for the lebotatoiy litem- 
/V pi« and the efficiency of Hj’ producflwnn 
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DLl and DL2 were tep irately detarmined. 
During ptobt descent, the ’He/ 4 He ratio 
waa determined from 3 -emu and 4-amu data 
obtained with die noble pi cell (NGC) 
temple (6). Hydrogen waa effectively abaent 
from dll* gai aample. The reiult, 1,1 x 10 -4 
(Table 1), Indicate! that He tn the jovian 
Mmoephcrc ha. a .mallet proportion of ’He 
then tha prototolar H. found in mataorlm, 
where J H./»H* - (1,5 t 0.3) X lO -4 (8). 
Evidently, the "protoeolar" He Itotopci in 
meteonm do not, in fact, reflect the values 
that existed In the solar nebula. 

With the He mixing ratio and HZ' pro- 
duction efflcienctea in hand, the HD/H* 
ratio wu detarmined from the DLl and 
DL2 data at 3 amu and 2 »mu. During DL2 
sampling, the tUctran en.igy in the ion 
toutce was reduced at timet from 75 eV to 
25 eV or 15 eV. Few and no He km*, 
respectively, were produced under there 
conditions, and tha ratios of die 3 -amu rat* 
to dia 2-emu rate provided an upper limit to 
HD/H] (aaumlng no production) of 
(1.1 ± CJ) x 10 -4 . The avenge value of 
HD/Hj, m. a lured from inlet 2 data, wa. just 
at thi. limit, The ratio* maaeured in Itiletf 1 
and 2 were elmott pt.cU.lv th. tam*. The 
ll C /*’C Isotopic ratio wu deteimined with 
high precision in DLl data where die con- 
tribution of NHj at 17 emu wai negligible. 
The ratio i» exactly that of tha tolar values. 

Reuonably good agreement w» found 
for CH 4 mixing ratios tn DLl and DL2, 
within the calibration uncertainty of 
£20%, Tlilt determination Indicate* diet 
atmospheric ”C U mot. abundant by a 
factor of 3 relative rn H, than in (olar 
lystem material. The maasurad mixing ratio 
of H,S In die well-mixed atmoaphere be- 
tween (land 1 1 bin Increased by a factor of 
4 with depth. There wea a tentative detec- 
tion of HC1. Water vapor and NH> 
readily mburbnl on metal .utfocei, and 
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their daniitlei in the Ion *ource may lag 
behind thcae In the gai itream. The me>- 
(urementt made while DL2 wai open thould 
be die moit reliable, but they were made 
juac after enrichment cell 1 (BCD wu 
closed. Some of the NHj and HjO detected 
in DL2 may therefore have been contribut- 
ed by outgauing of thi* odwebed gu tn th. 
pretence of abundant Hfend He (P). Hence, 
for rha pruent, only dppre limit* can be Mt 
for these molecule*. The constraint for 
KH,, 16 time* the tolar value, it weaker 
than the exlating estimate, 1 to 1 J ttmu 
the ioiar value, from analytit of Jupitar’i 
microwave emlMlon ipacoum (10). The 
limit on water vapor, 0.2 time, the tolar 
value, i. ngn if leant in vi»w of the tolar or 
greater value* found for the other volatile* 
C, N, and S. 

Two- and three -hydrocarbon spade* ap- 
peared to be present. The direct leak* 
showed no sign of any heavier hydrocar- 
bon*, Upper limits for the mixing ratloe of 
three heavy hydrocarbon specie* are ~1 
ppb. 

Our H./H] ratio (0.156) agree* almost 
exactly with the ratio of 0.137 ohtalnod by 
the Helium abundance detector (HAD) an 
the Galileo probe (ID: tome of the tmpli* 
cation* of thi* measurement are dltcuned 
elsewhere (II). This value I* considerably 
larger then the value of 0.11 obtained from 
Voyager data (12), but is llttla (if any) lam 
than die value for the present sun. Tire 
pru mauler value, dedu c ed from evolutionary 
modal* of th. tun, i* probably a more rele- 
vant itandard and t* 1854 grantor (11). It ii 
widely believed that Saturn 1 ! atmosphere, 
with a He/H] ratio about one-fifth of the 
prototolar value, hu been depleted of He by 
rein-out of He droplets in the interior (l 3, 
4). The OPMS end HAD results eg™ 
l i i l i iln V i [ I I r 1 >n auggasting that 
the same process must operate at Jupitar, 
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though not to the Mine extent Thla Ion of 
He is * likely explanation of the Ne deple- 
tion observed by the GPMS, if the solubility 
of Me In He is m large m has been suggested 

(15). 

The value of (i ± 2) X 10~ 5 deduced 
for D/H from HD/Hi is hither thin the 
ratio* derived from Tsarth-based spectro- 
scopic obiervaeiom (J6, 17)i (2 £ 1) X 
10 5 from CHjP/CH 4 and(2il)x 10“* 
from HD/Hj. It * grata well with the value 
found in a recent analysis of Voyager infra- 
red ipeetra (l$h (3-8 ± 0.5) X 10" The 
jovien D/H ratio should be the ww 9% the 
protnmlir or solar-nebulit value because 
the nebula was the source of Jupiter’s hy- 
drogen, [The terrestrial ratio is 1.6 X 10 .] 
The sun converted ell Its D to J Ht early in 
its history; thus, Uciss (B) argued that the 
original D/H in the solar nebula should be 
equal co the difference betwmi die 'He/ 

4 He ratio* in the solar wind and in meteor- 
itfii which should contain nebular He (19). 
Gain's value for pmmanlar D/H Is (2.6 ± 
1.0) X 10’ 1 . and this overlaps the GPMS 
valuer Substituting the new determination 
of protoeoiar ’He^He - (1.1 ± 0.2) X 
10 -4 for thi miceorttic one. the same argu- 
ment would yield D/H ■ (3 s 1) x 10“\ 
a value that is closer co that derived directly 
from HD. It ii wVser drat D/H on Jupiter l* 
distinctly hightr than the value of (1.6 ± 
0.1) X 10"^ found In local Interstellar H 
(20), as expected from the steady destruc- 
tion of D in the galaxy by stellar nucleosyn- 
thesis during the last 4.5 X 10* years. 

Earth-based imaging (21) snd the re- 
markably small number of cloud particles 
detected by the nephe lo meter experiment 
(22) indicate that the probe may have de- 
scended into an unusually clear part of the 
jgrvian atmosphere. An obvious explanation 
U that the region was one of subsidence, 
like most clear regions rm Earth* Our H : 0 
abundance value, then, may be Iras than 
die planerary average, although this result is 
unlikely. In a cloud, the relative humidity 
should be just over 100%, and this is con- 
sistent with die comprehensive analysis of 
Voyager spectra by Carlaon it al. (IS). The 
corresponding profile of water vapor has » 
scale height of about 3 km* If for example, 
the relative humidity U reduced by a down- 
draft from /D0% to 90%, the amount of 
subsideflafe required is 3 ln(90/100) km. or a 
men r0.3 km. Our measurements were 
made at for greater depths than the expect- 
ed condensation level (5 or 6 bars); to 
perturb the humidity at such a deep level 
would require a downdraft extending over a 
range of well over e scale height, for more 
than needed simply to clear out the clouds. 
There remains the possibility of a global- 
scale circulation, perhaps upward at high 
latitudes and downward it low. 

The large atmospheric abundance of C 


suggest* that sources other than the toe of 
the solar nebula contributed a significant 
share of the volatiles in the envelope of 
Jupiter. These elements and water should 
have been in the condensed phase in the 
neighborhood of the accreting Jupiter One 
possible explanation for a deficiency of wa- 
ter was proposed by Gautier and Owen (j, 
4); most of the plane tesiroali accreted to 
die rocky, icy core of the embryo planet 
rather than being collected after die planet 
began to amass Its gaseous envelope from 
die neighboring soler nebula. Water vapor, 
being the least volatile, remained close to 
the core, whereas gases like CO, CH 4 , and 
N, could mix with the ncbula<derived en- 
velope, enriching it In C end N. Although 
this scenario could reproduce the observed 
abundances, it must be tested by rigorous 
modeling. If then was a luge mass of late- 
accreting planets rima Is, the low relative 
abundance of atmospheric H z O requires 
them to have been rich in caihrmacentM 
materials and deficient In O. It has been 
customary to assume that outer solar-system 
planeteiimsli resemble the le*st»altered 
carbonaceous chondrites that fall to Earth, 
which contain more O than C (23), but 
there are many that an much drier. The 
mean d trial ties observed for many satellites 
in the outer soler system have been ex- 
plained by assuming e rock: Ice ratio of 
order unity, but perhaps much of the low- 
density material is carbonaceous rather 
then icy. 

If our results are representative of the 
entire planet (24), they preclude models tn 
which several Barth masses of late-accret- 
ing, water-rich planeteslmals deliver vola- 
tiles to the juvian envelope ( I ). There U no 
sign of the dense witer cloud, or the corre- 
spond mg vapor, invoiced to explain waves 
spreading from the Shoemaker-Levy 9 im- 
pacts (25). Water vapot observed In the 
Impact plume* probably came from the Im- 
pacted, not Jupiter, as suggested (26). 

mfimncmanp noth 

i . o. J, ateransen. Amu. aw. ton*. M. 10. 
25? 11012); J. a. tote* aid P, towMmar, in 
Origin sndftouwi wasnetey ** Samara aido- 
I fifm m , *• K. Afrqe, J, 6 tote*, M 0 MmNmi, 
Kell. (UW. of After* Press, Tuoson, 1080), p. 564; 
WBHubbero&d.O, 980. 

1. J, I.LunfneandO.J. 6tevontOT,AtNpfy*. / Buefr 
SI 403 (1 068), T. Owen and A. to -Nun, taarue 1 1 0, 
216 flOWk 

а, O. taitte anU T, Ow, m Ortgh end fwtattan of 
Ptoeey toSMtoAerraeoherei. I. K Atreya. J, 
B Pollack, M, 6. Matthews. ids (Urfe of Artone 
tows, Tvceon 1*60), p. 46?, 

4, K. 0. Hd end H. P.Lsraon, towns <6. 101(1001). 

б. H.».Msm*irefef.,*oeeticiAte,6ail1(iS02). 

6 , each Wat dso Ms gee senate loin enrairroent 

shembor, in whtoh o parous eertran meteria eb- 
s Dib e dqo m pte ftytiooerbcnsmdhei ^newsgas- 
er The noon oases not adsorbed by act were 
captured and M te the ***** n on N0C tom 
wfson H* was eikhnatid by a getter pump, tni 
pfooedirtlwraradtratMehMotdetictsbityt^^ 


teeter ef 10 (te 0.1 pal par NMer by volume}. The 

oeeee not absorbed SiKIvraraaiH pumped My, 
end toe adsorbed am wvt later ratoed by neat* 
hg era tel and endyto These meaeunmento, to 
two mearofemsrta of teBkground wve eerrled out 
iter Inlet 1 vras otand and (tore Hti 2 was 
opened, IC2 wee M 6em Wet 2 when the pres- 
sure wee M to 0 brae. The two tochmer t oeU 

tea a te eAebte te m m m m m •*- — -4 

raeeomeeenemty wpervorooerooneenome 
heater neoie gee* to ee much m e totor of 800. 
The svtiem wee oonetod by e reed-on*/ memory, 
incremented each 19 a tough 8102 itape. Many 
scene were nrteOe m Nqr* m e ese e through the 
rinoe of 2 to 160 aiu reotirtno 75 •: other ee* 
chances coeradcrSyatiM id mrases.Ttetywn- 
tomge woe 10*, 

7 I. Anders end M Ommm, flenrten ft n smarte n 
Act* Ml 18T(1W» Thm itunAnow, oornmonly 
«vtm to at "Mtr' « '‘icNr^yiwn'’ ummmo> 
ee, ere based on a oeverador of eNtdonoo Pom (he 
sun ttsdf to tom prtnMue meteortes* 

a J. Gads, »***/*» of *>« fldmwxe, N. 
Prantaee, K. Venpiwd-tem, M, Cued, Bde, (Com- 
brtdoe Unlv. Prase. Oembrdpe. t W8), p. 66. 

0. The enridimenk oeSe were heeled end the evoked 
g eeee meeeured be tw een 6,0 end 7 bar and 11.0 
end 14^ bar. 

10. A Merton, R Courte D- Oeurar, A. uoombe, ios- 
rur 41, 410C 1 660k l>te Peter andfoT. SAute, fetal 
01 143 (1086). 

11. II. von Zihn and D M. Hunan, Safmot so, m 
(I860), 

12. 0 J. Conrath. D. deter, R A Hend. J. 8. Hnm. 
stall, Atrophy Ji m, 00? {1064], 

13. aarauoncMM./eMurarioei (w). 

14. o* Ote ve n ae n end B. I. tefodter , Atertehya. J. 

$up0, 26,221 (lOTTfcted, p.230. 

1 5. M. 8, Poueton end & J. St e raneon, ID! ikprrig 
eun *>76, 040{i wag 

16. 0, deter end T, Oran, It 0), p.004, 

17. w. k amifh, w. v. achempo. k. h. tee, Aem> 
ptry*J, 10, e6V(1f*R 

ia. a f. Coriflon. A. A. Lade, w, B. Poaaow. JL Qeo- 
phyw to e, IS, 6861 (1603). The D/H ratio Quoted In 
me riport nee been eUiueteU tu itect egnUy 
leroer moomement d CH 4 abundance. 

10, The soier 'Mndveiue Is (4.6 2 0.4) x 10" 4 , end the 
meeontconeeOAt03|x ic~ 4 . 

20. J.LLbiteyteeA.Ateephy*, J. 401.004(1 000}: eee 
esc X, Uniky teat, k\Uf* Sknmt Abunawioaa, P, 
Crane, Ed. fo^tioer-tortsg Bern, 1099. p. 21 6. 

21. O. CHon or 4 wm.ww (1000, 

22. B. Rag enl. D. I. Coton. P, Awn, K. A, tooes,ldd. , 

p. XXX. 

23. H, B, Wnk, OuMten Coemoechim, Acte 6, 270 
(1069: a. Meson, Stood— Oteey. Hew York. 
1069. 

X4. A pro knwpitecn Of the obeerved verwUcna In 
abundenoee coutd toterdeduely ante dbtoncee 
h lh. wU9My q( mMh In mMMe H. Afitnd- 
m««y or Ji4Mr> h • k m mauac pn«M. 
Spr inateno*. S the re * teg rede of a minor oohetiuara 

m me mtieotfor oetaftea plane « rrocb drger then 
ire solar value, tree m msena raao e me motto 
toHorraay bo eighSQeieSer than tie eater vtiueor 
vtoevaraa. 

25. A, P, hgeriOP end K Konomorl, Atafura 874, 709 
(1069 

26. A.L.8oraoueeel,toui < lrprasfe 

27. We her* J.Cooiay tar We etertseilnetor onlmin- 
ever, R, Lull end t. Tyte for the mechanical dtegn 
oontributione, A. Ooet tv Ne hdp hihe anrtohmerd 
odi deep, ft. AM, Hi Powen, mt H. Mends tor 

the predeton —wi tej wetatog endjtetfSnin g , end 
R Araey end H. B ed teilar thi dete f ontoi a e aembhr 
eti testing at Codded Bpaoa Plgra Center, The 
oontrtouiione ta thetonord ■ eyeusn Oeton ci i. 
Stock, j Certetil. J. tor, J. Merer, end W Phkus 
a the Untordty of Mtogen era grtteuly edkwM- 
edged* T *» W teto ctote were torfcteeU 
et tto totertH FtenMe Aerm tpaee QMtidn it Veley 
Forge, PA snd the wtew te ee tor me g — a empUng 
sytram weraoeetofl rod totoed by am m» 
toee in Otonh. CA. Ite getter mstoi 

wee prodded by the MB Oerter ot Miien. Italy. M. 
wong or the unherray or Msysgen p erttop ra ed in 


SCdNCB • VOL. eee ■ eee 1006 


3 




04-24- 199G 1G:0H 


30i HOG 17GC 


NAG A CCrC CODE 91G 


P.23 



tht dw anH/M. wi mo think th# Omo * obt 

Profact 4 NASA Arrm A iin rcr Ctntir 

md wa oftoiirV ■e Anada d p i tha oor ti touttona of 


A. VVtwW, C. Soboefc, md fc Mik* to th«ir uttarti 
during tha davdopmant, Ipi OKra ll Hfgr*len, and 


tilting phwM IS Stooh 1WS; aco«o<$d 17 /^rfl 
IMS 


d 


cnrvrr * V<*M • +++ 1 r ' ' 




07-09- 199C 0C s 14 


301 20G 1 73 G 


NAG A GGrc CODE 913 


P. lO 


CoS ? Afc '/?<• 

Chemical Composition Measurements of the Atmosphere of 
Jupiter with the Galileo Probe Mass Spectrometer 


H. B. Niemann', S. K. Atrtyi\ G. R. Ctrignan 1 , T. M. Donahue 1 , J. A. Haberman'. D. N. 
Harpold 1 , R. E. Hartie 1 , D, M, Hunten\ W. T. Kasprzak’, P. R. Mihaffy 1 , T. C. Owen 4 , N. 
W, Spencer 1 


1. Goddard Space Flight Center, Greenbelt, MD 2077 J USA 

2. Department of Atmospheric, Oceanic and Space Sciences, University of Michigan, 
24S5 Hayward Street, Ann Arbor, MI 48109 USA 

3 . Lunar and Planetary Laboratory, University of Arizona, Tucson , AZ 85721 USA 

4 . University of Hawaii, Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, HI 
96822 USA 

5. 12013 Remington Drive, Silver Spring, MD 20902 USA 


ABSTRACT 

The Galileo Probe entered the atmosphere of Jupiter on December 7, 1995. 

Measurement* of the chemical and isotopic composition of the Jovian atmosphere were 
obtained by the mail spectrometer during the descent over the 0.5 to 20 bar pressure 
region over a time period of approximately 1 hour. The sampling was either of 
atmospheric gases directly introduced into the ion source of die mass spectrometer 
through capillary leaks or of pi which had been chemically proceated to enhance the 
sensitivity of the measurement to trace species or noble gasee. The analysis of this data 
set continues to be refined based on supporting laboratory studies on an engineering unit. 
The mixing ratios of the major constituents of the atmosphere hydrogen and helium have 
bean determined as well as mixing ratios or upper limits for several less abundant species 
including: methane, water, ammonia, ethane, ethylene, propane, hydrogen sulfide, neon, 
argon, krypton, and xenon. Analysis also suggests the presence of trace levels of other 3 
and 4 carbon hydrocarbons, of propane nitrile, phosphine, hydrogen chloride, and of 
benzene. The data set also allows upper limits to be set for many species of interest which 
were not detected. Isotope ratios were measured for *He/ 4 He. D/H. '‘C/ 1J C, "Ne/”Ne, 
*Ai/ w Ar and for isotopes of both Kx and Xe. 
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INTRODUCTION 


The determination of Ihe composition of Jupiter’s atmosphere and the comparison with 
solar abundances is of interest from the point of view of undemanding both the nature of 
the formation of this planet through mechanisms such as direct capture of nebular gas or 
infiision of icy planeteilmali and the subsequent evolution of Its atmosphere by 
precipitation of constituents Into the interior. Although remote sensing from Voyager and 
other spacecraft and from ground based observations had coos trained the relative 
abundance to hydrogen of helium, methane, water, ammonia, phosphine, and other 
species prior to the Galileo Probe entry, sampling below the cloud layers was only 
possible for many species with the insitu measurements enabled by the Probe. Isotopic 
determinations which had been obtained by earth based spectroscopy for some species 
such as deuterium [Smith, et al., 1989, Carlson et al., 1993] could be measured for a 
broad range of species with an inaltu mail spectrometer measurement. 


Preliminary analysis of the GtUleo Probe Neutral Mass Spectrometer (GPMS) data 
[Niemann et al., 1996] gave isotope redos for the noble gases, D/H, and '’C/’C, and * 

abundance ratios to hydrogen for helium, methane, neon, argon, and hydrogen sulfide as X 

well as upper limits to the mixing ratio for water, ammonia, krypton, and xenon in the 9 <y^ 

to 1 2 bar region of the atmosphere. £cntinu>nt analysis of the flFMS data since this 
preliminary report has established the presence of the additional species e thane rethyteriET^V^tf , <> ^ 
and propane, suggested the presence of jaed rog e n cM ortde, benzene, and propane nitrile, 
and demonstrated an unexpected and substantial increase in the abundance of water, 
hydrogen sulfide, ethane and other specks with increasing depth into the atmosphere over 
the 1 1 bar to 23 bar pressure regime. It has been suggested treys ??, 1 996s Owen et el., <>Ji ' 

1996] that this increase may be due to the probe entry into l 


large scale vertical atmospheric circulation structure which has an atypical abundance 
many species compared to the global avenge. This depletion it attributed to precipitation 
of gases moving upward through the cold region of die adjacent member of the 
circulation cell. 


EXPERIMENT DESCRIPTION 


<sl.fi 


Instrument Description 


The instrument has been described [Niemann et al., 1992]. The instrument consists of • 
quidrupole analyser mass spectrometer with a mass range 2 to 150 amu. The ion source 
had a magnetically confined electron beam emitted from a hot filament. The electron 
energy was nominally 25 eV with occasional 25 eV and 15 tV measurements. The ion 
detector was an electron multiplier operated in a pulse counting mode. The ion source and 
analyzer volumes were continuously pumped by a miniature sputter ion pump and by 
chemical getters. The gas sampling system allowed gas horn one of two direct leaks 
(DL1 or DL2) or from enrichment cells (BC1 or EC2) to be introduced into the ionization 
region of the miss spectrometer. The operation of one of the enrichment cells was 
optimized for the detection of rare gases (RG). Several background spectra were taken. 
The gas sampling system and other instrument elements are shown schematically In 
Figure 1 . Figure 2 illustrates the time sequencing of the direct atmospheric 
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Figure 1. The OPMS inlet lines, enrichment cells, capillary leaks, micro valves end other instrument 
elements tit schematically illustrated. From (Nlemmn ct ■!., 1992). 


measurements, the first 
background measurement, and the 
sampling of the atmosphere by the 
enrichment cells over the pressure 
and temperature ranges 
encountered by the probe. 


The gas sampling system allowed 
atmospheric gas flowing through 
erne of two inlet lines to be 
sampled directly or after chemical 
processing by one of two 
enrichment cells. The later 
measurement! were designed to 
provide a cleaner noble gas 



Figure 2. The sequencing of the 2 separate inlcti and the 
enrichment cell sampling times are illustrated. 


Temperature (K) 
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measurement end in increased 
signal to noise level for tnce 
species. In addition background 
measurements were taken at 3 
intervals during the descent 
sequence. The measurement 
sequence over the 6863 
integration steps during which 
data was transmitted and the 
estimated effect on the pressure 
in the mass spectrometer Ion 
source is illustrated in Figure 3. 
The detailed step number and 
the time at which each 
sequence was initiated is 
summarized in Table 1 together 
with the approximate temperature 
and pressure at the start of the 
sequence as derived from the 
Atmospheric Structure Experiment 
[Sief, et al., 1996]. 



10 * 

itr* 

i<r r 

ic B 

nr* 


Figure 3. The descent sequence is illustrated by 
showing the approximate ion source pressure at 
various periods during the descent. Background 
measurement period! and enrichment cell periods are 
also illustrated 
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all source inlet valves dosed 

RG (rare gas) 

gPl 

1092 

268 

4.6 

V2.V3 closed, V4open 

BC1 (enrichment cell 1) 

G 

1238 

281 

5.4 

VI , V2, V5 closed, H3 on, V3, 
V4open 

B2 (background 2) 

3097 

1561 

306 

7.3 

all source inlet valves closed 

Dl 2a (direct leak 2a) 

3482 

1753 

321 

8.6 

inlet 2 open, all other source 
inlet valves closed 

BC2 (enrichment cell 2) 

4442 

2233 

333 

12.1 

V8, V12 open, V7, V9 dosed, 
C2 heated by H6, signal 
through L4 is added to that 
from L2 which has no shut off 
valve 

DL2b (direct leak 2b) 

5331 

2677 

381 

15.7 

V8 closed ' 
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Initfumant Performance 

Although leveril of the probe experiments were affected by the greater than expected 
thermal coupling to the probe interior of the Jovian atmosphere fYoung et al.. 1996], the 
data set from the mass spectrometer was relatively little Impacted. This wit due, in part, 
to the protection of thermal inertia of the internal pressurized housing which had been 
implemented to avoid electrical discharge at certain pressure levels during the descent 
Temperatures in the Instrument electronics varied from approximately -IOC severe! 
minutes into the descent to approximately 70 C near the end of the probe mission. 

The response of the detector to changes in the sampling system for several selected 
masses over the course of the descent is illustrated in Figure 2. The signal levels during 
tha background measurement period! and the enrichment cell periods indicated that 
valves and heaters performed as expected 


1140 

1«+7 

i'+f 
le*5 
1e+4 

1*40 
1*42 
1#4l 
1a+0 

0 1800 2000 3000 4000 MOO 8000 7000 

step number 



PI jure 2. Detector counts for selected mass values are shown over the entire descent 
sequence. Relatively large toatases in the count rate at masses 18 and 30 after EC2 are 
evident The counts have been corrected for detector dead time effects. These correction* are 
valid over the entire descent for the counts shown at 30 emu, 18 emu. and 16 emu. However, 
in the later portion of the DL2a region, the EC2 region and the DL2b region the detector is 
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nearly saturated at 2 and * amu and for these measurements the form of ihi« corceciioo it not 
valid. 

The variation in ratios of the countt/period in the DL1 and the DL2 region for the 
different males illuttnte the need for independent calibration* for the two capillary leak*. 
This it true, in general, for each leparate species. DL2 introduce* the gas into the 
ionization region in a direct beam path from the end of the capillary leak and a portion of 
this gai la ionized before it tharmalizaa with the wells of the chamber. The peth to the 
ionization region from DL1 is more indirect because of the requirement to dose valves 
viOa and VlOb after the flm segment of the measurement is complete. 

COMPOSITIONAL ANALYSIS 

Summary of previous Rrenlti and New Findings 

We have previously reported abundances of 0.156 ± 0.006 for ‘He, (2.1 + 0. 1 5) x 10 : * for 
CH 4 , (7.7 +■ 0.5) x 10 4 for H,S and Upper limits of (3.7 + Q.15)Kia* for H,0, (3.5 +■ 0.3) x 
Hr 1 for NH P (8J +0.4) x 10 4 for Kr, and (5+2.5) x 1 O’ for Xe. The helium to hydrogen 
ratio it in excellent agreement with the value derived from another Probe instrument the 
Helium Abundance Detector (Von Zahn, at al., 1996]. Methane abundances were derived 
from measurements in both the DL1 and DL2* region. Hie minimal Interference at mass 
17 in the DLl region allow the ”C/*C ratio to be established at 0.011 [Niemann et al., 

1 996], identical to the solar value. 

The upper limit of 3.7 x 1 O'* for water previously repotted [Niemann, et al., 1996] was 
derived from the region just below 1 2 bar (prior to step 4442 where the enrichment cell 
measurement was initiated). This value represents a limit of 0.2 times the solar value. 
However, the data shown in Figure 2 shows that the water abundance increased by the 
time the enrichment cell experiment EC2 had been completed and direct sampling of the 
atmosphere resumed at 15.6 bar. This rapid Increase In mixing ratio is also observed for 
several other species and is consistent with the view that atmospheric circulation effects 
created a large region substantially depleted in condensable species at the higher altitude 
regions of the probe entry she. Another example is the variation In mass 30 counts whh 
step number, also shown in Figure 2. 

The DL2b sequence provides data from the deepest pert of the atmosphere sampled 
before the instrument began to degrade due to the extreme conditions. High quality data 
was obtained during the DL2b sequence overthe pressure range from 15.6 bar to 
approximately 2 1 bar. This region shows the greatest number of chemical species in the 
spectre. In addition to the species listed above we have identified ethane, ethylene, and 
propane in the spectre. Our preliminary analysis slso suggests the presence of trace levels 
of other 3 end 4 carbon hydrocarbons, of propane nitrile, phosphine, hydrogen chloride, 
and of benzene. 

2m EasauUL Laboratory Calibration 

The calibration of the flight unit in 1985 sampled a variety of chemical species postulated 
to be present in Jupiter’s atmosphere. Several of the species identified in the GPMS data 
were not included in these calibration experiments and we are presently continuing these 
experiments with • spore GPMS unit. This calibration Is necessary for quantification of 
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mixing ratios since the sensitivity of the Instrument to different gases varies considerably 
from species to species due to the gas inlet and pumping systems. 
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The Galileo Probe Maas Spectrometer: 
Composition of Jupiter’s Atmosphere 

Huso B. Niemann, Sushil K. Atraya, Qaome R. Carignan, 
Thomas M. Donahue, John A. Haberman, Dan N. Harpold, 
Richard E. Hartle, Donald M. Hunten, Wayne T. Kasprok, 
Paul R. Mahaffy. Tobias C. Owen. Nelson W. Spencer, 
Stanley H. Way 

The oompo*Won of tha Jovian itmoephere from an altitude merited by an atmoaphartc 
preeeure of 05 to 21 barawaadatarmlnadbyaqLiadrupoiamaaiapactromatarontha 
Oaato pwoa. Tha mudnfl ratio of Ha (heMum) to H, (hydrogen), 0.1 68, ta float to ore aoiar 
ratio. Tha aounoanoaa of mathana. watar. argon. noon, and dihydrogen aUtda «•» 
maaaurad; krypton and xenon wara detected. Aa meaeured In tha Jovian atmoephere, tht 
amount of carbon la 8.9 tlmaa tha aoiar abundance ralatlvoto H B , tha amount of etlfor 
la groator than tho aoiar abundance, and tho amount of oxygon la muahlooa than theeolar 
abundanea. Tha naon abundenoe oomparod with that of hydrogen la aboid an ordar of 
magnitude leee than the aoiar abundance, laotooic ratio# of carbon and the treble gaaa* 

art conalatant with aoiar vggaa. Tha maaaurad ratio of dautartum to hydrogan (pAS of 
(9 + 2) x 10"* Indteataa that thla raSlo It greater In aotar-ayatam hydrogen Start In local 
Interetelar hydrogen, and tht *HarHe ratio of (1 1 1 0.2) x 10- provKSoa a now valua 
for protoaoiar (aoiar ntowaj hoi lum leotopca. Together, wa EVH ana •mf 4 ** ratio* am 
conatatart with oonvaralon In tha avn of protoedar dautartum to preaent-day *Ho. 


Determination of the comportnot* nf Jtirrf- 
tar’i atmosphere ihould oonttrain the rela- 
tive Importance of direct contribution* to 
tha atmoaphere from the tolar nebula itielf, 
on the one hand, ml from large Icy or 
jocky object! pteient in the early outer »Uf 
•ytcem. on the other. The contribution of 
(hew “pkawtetlmala* could have been In 
foe form of an early etmoaphem around 
Juptcer’i primitive lolld com or In the form 
of volatile* they carriad to the planet after it 
befan to acquire ica atmotphtta of H and 
He from the tolar nebula (1,2). The degree 
of reeemblance between die nm't atmo- 
qfotre and that of limiter In the abundanc- 
es relative to H or elementi tuch a* C, N, 
and O ihould be decMve in provldlnt the 
M|ulnd concertina. 

Before the direct. In tint maaauremsntt 
reported hen, ntmu •peeaoeeople tenting 
from Berth and fare i p ec ecw ft had Indkar- 
ed that (he telathra afaundancM of C (in 
CH«) and N (In NH.) were graatar than 
die win abundances (3, 4). Fwaphlne and 
water vapor had bten detected, but thetf 
abundan ce* below the doudi were uncer- 
tain, and HtS and die noble gaeci odier 
thin Hs Km net bocr dtfittltd tt til* 

K I. Hmmn J. a. Hfesrrwi. D. N. HsrpoW. R-C 
Hint, W. T. *Mpn«< P. n. MMy. Ok *w*y Job- 
dvd 8 c** Crnm. (MS, MO S 07 T 1 . UIA. 
I K. *r *8 0 . 1 * Carter, T f M. Drahut, Un*sr*ytf 
MOSQW 24M HpyMrt Avwfc Am ArUP. Mi *018* 

0. M Hutu* UrtusnSy Awn* Toetcn, AIWT21, 

HonoUu, H UI A. _ _ 
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determining the itotopic redos of their con- 
stituent elemeoa (5). Signal* ware recorded 
from over 6000 individual valua* of man. 
co-charge ratio (m/t). The data returned by 
the tnetrument have led to the diacovety of 
tlx are atinuapltcifc umreitueao ad (hs 
acwuRiMitt of numcmui abundance* and 
isotopic moot. The gat templing eyrtstn 
admitted rhe jnvian irwwylir w through 
low-conductance leak* to the Ion tource of 
a quadripole matt analywr. The ion tource 
wet pumped by t getter and the analyter 
volume pumped to a much lower praam 
by a getter backed by i eputtar-ion pump. 
One inlet (direct leak 1: DL1) waa open 
6om 0.52 to 3.78 ben: the other Inlet 
(DU) fmetioned from 8.21 ban to the end 
at about 21 bon (<)• 

Rfun 1 dutwt a ample maa qwcmim 
from tha 3-bar region at 228 K. before my 
significant watar cloud wa* ejqwcad, end 
another one from die 11-bar region, where 
the 330 K temperature unite* that the 
probe wa* (hr below any water eondentatton 
lrveU Curopwiaun uf tire twu awetre *how* 

a lane increase in die w end direct leak ef 

HP and NH,. In firot, the design «f <ho 
Irak 1 templing tyierm. with • relatively 
finlttve dita to replace there tanalMtig Iona vacuum path to tha Ion tourca.p«- 
chiet to the tourcet of volatile compoundi -• dude* a retuldye detection of the tunace 
In Jupltert^rmoepherc were noteyailaH^mctive ipedet HP end NH. in the iM» 
before tW*3eiceni of A<^*ar^fij2_.*phertc region templed by mu leak. Raw 
An acowate direct meuuronent of the count! thown In mete tpectre mutt be core 
abundance! of He and the other noWe garea recred for the counting tyttam deedclme. 
u Important for undemanding how Jupiter The cfTicienda t* “ re loo source end are 
hai procured to conatitucn» Abutvienec* ^/pumping ry*tcm vary from qreatre to *pe- 
lere than dreaa af dhw tolar valuei would/ clea, and dre w 


ilgnal the formation, in the metallic Interi- 
or of tht ptonat. of a repartee He phaae. 
tome of which ha* precipitatad deeper into 
the Interior. Remote renting had indicated 
a modett depletion of He, but AS indirect 
reault la much leat mtufoctory than the 
direct renting by two of tha Oalilao innu- 
mcna repotted here. 

Tk* OmKUo probe mis spoctratuto 
(OPMS) waa daatgnad to maa re rr a foe mi*, 
mg rati oi of major and minor qwclaa while 


Mtwltivity freon, afoiA eta 

i t ffr w nt for the two direct teaki. mutt be 
todtvidually determined. In moat ceaae. we 
reed preflight calibration*, but for Ac moat 
layortant tpecla* we were able to do mow 
recalibretioni on tha (pare Inrenanant. 
Wtth the direct kaka, the background level 
wre typically between 0 and 3 count*. Two 
background mcaauremcna damonmated 
that foe Iniaument strived at Jupicat in a 
vory clean itae* and quickly lore in memory 
d almost all pare* contributed by CL1. 


Tame 1. Mxfig retire io H, and woioplcridoee>fw £4re> mi: wre»tmu. Ban* epeceammereaw*. 
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TabU 1 show* tha m bring retire (num- 
her densities relative to H,) of weeks that 
were detected and for which tellable mea- 
surements could be made. Error estimate* 
are toughW lfr and represent the dan seat- 
ter about the average value. Some of the 
obeerved mixing retiae can be compered 
with thaec for the sular *y»»*m aa compiled 
by Anden end Grcvcssc {7), commonly 
referred te at the telex ebundmeet. Alto 
shown ere result* of earlier tnalytef bared 
on Voyager or Barth-baaed spectre. 

The jovian H*/H> ratio. 0.156 £ C.006, 
la the a rere ge of 14 Individual tneaaure- 
manta from DLl (0.137 t 0X06) and 4 
from DU (0.131 £ 0.004). All the other 
noble gaiea were detactad. and they ere 
ckmcMially fractionated. Tbt N« ratio to 
it about an eider of magnitude lew than the 
io lor ratio, wKamat far Ar the ratio It slight, 
ly pee ret. Kr and X* ware clearly de tected! 
upper limits on chamfering ratio* an rttpac- 
tively 5 and 50 timee the eolar value*. For 
all four of the heavier noble reata. me 
iaotopic ratio* we cloee to die tolar value*. 

The ‘HtfHe and dia D/H nfloe were 
determined from an analysis the 3 emu 
(auaiiic maw unit) data. This analysts It 
complicated by the need to correct far H/ 
produced In the ion mme end to disrin- 

J tiith between die contribution! of HD and 
He. H, 4 U generated in one of two wayti 
by dissociative I on lotion of CH 4 or by a 
two-body ion-molecule reactor between 
Hj + end H,. The efficiency of H, + produc- 
tion could DC determined during the second 
background study when little Hr end He. 
but a considerable amount ef CH,, aaiitad 
in the ion tnnrce and CH. wu the only 
important source of 5-emu ions. During pre- 
flight calibration, when the Hapressure in 
the GPMS WM varied om i < 4rge tenge, 
the HD/H, retfo for the kbotatoty HtaO- 
1 pie and the efficiency of H,* producflBrfln 


DLl end DL2 were separately determined. 
During probe descent, the *He/*H* mdo 
was determined from 3-tmu and 4 -emu data 
obtained with the noble pi cell (NOC) 
sample (6). Hydrogen waa effectively abient 
from this |u sample. The result, 1.1 x 10 -4 
(Table 1), Indicate* that Ha in the jovian 
Mmoaphcrc ha* a amaller proportion of *Hc 
than m* prototolar He found in mareoriere, 
where J H./*Hs - (1.5 * 0.3) X lO" 4 (*). 
Evidently, the "protoeolar” He isotope* In 
meteor res do not, in fact, reflect the value* 
that existed In the solar nebuk. 

With die He mixing re no ind H,* pro- 
duction efftctencie* in hand, the HD/H, 
ratio wu daremined from the DLl real 
DL2 data at 3 anui and 2 amu. During DU 
mapling, the alaatson anergy In the ion 
reutea wa* reduced et time* from 75 eV to 
25 eV or 15 «V. Few end no He Iona, 
respectively, were produced under there 
condition!, and the retire of die 3-amu tern 
to the 2-amu ate provided an upper limit so 
HD/H, (amumlnc no H, 4 production) of 
(U ± 0J) X 10 -4 . The avenp value c# 
HD/H,, measured from Inlet 2 data, w» jure 
at this limit. The ratio* maaeurad In Inlet* 1 
and 2 were almost precisely die same. Tha 
lJ C/°C isotopic ratio wet determined with 


high precision in DLl data where the con- 
tribution of NH, at 17 emu we* neglltfbl*. 
The ratio is exactly that of the solar value*. 

Reasonably pod agreement wa* found 
for CH. mixing ratio* in DLl end DL2, 
within the calibration uncertainty of 
£2014. Till* determination indicate* that 
atmospheric 1J C It more abundant by a 
factor of 3 relative tn H, then in solar 
tyetnm materiel. The measured mixing redo 
of H,S In the well-mixed atmrephere be- 
tween 8 and 11 ban Increased by a factor of 
4 with depth. There wee a tentative detec- 
tion of HCL Water vapor and NH, 
readily aihorbsd on metal surface*, and 
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their daniitiet in the km aouroe are kg 
behind theta In die gas stream. The mea- 
surements made while DU wa* open diould 
be dm most reliable, but they *« “*de 
just rfter enrichment ceO 1 (BCl) waa 
cktad. Some of the NH, and HjO dneted 
in DL2 may therefore have been contrib* 
cd by outgasting of tbit adsorbed gat In the 
pretence of abundant Hfod He (9). Hence, 
(nr tha present, only dppre limits on be let 
for there molecule*. The constraint fa* 
NH„ 16 timai the solar value, b weaker 
than the existing tatimata, 1 re I J tom* 
dm soiar value, from remlyai* of JagsrerT 
microwave emission spectrum (J0L Tha 
limit on water vapor, 0-2 time* tbs solar 
salve, i* significant in view of tbs rek* or 
gre a wr values found for the other volatiles 
CRandS. 

Two- and three-hydrocarbon spades ro- 
peered to be present. The direct leaks 
showed no sign of any heavier hydroeer* 
bone. Upper limla for the mixtag ados of 
these heavy hydrocarbon spectre an -1 

ppb. 

Oik Hs/H, redo (0156) agree* abuse 
caaady with m ratio of 0.157 obretood by 
cho hstiian abundance detector (HAD) on 
the Galileo probe ()l)t soma of dm impli- 
cation* of this measurement are dbcusred 
elsewhere (11). This value b considerably 
larger chan dm value of 0.11 obminml from 
Voyager dam (12), but u licda (if any) km 
than tha value far the prarent suss. The 
prouaular valus, deduced from evolutionary 
models of dm sun, is probably a more sek- 
vant standard and it IF* graatre (ID. It b 

widely believed that Saturn's atmorehere, 
with a He/H, ratio about ooa-ftfth of dm 
procosolsr value, ha* bean depleted of Ha by 
nln-out of He dropkta in dm uirenor (13, 
14). Tha OPUS and HAD luaula agm 

***Saisb >li* Vwysgcnfow in suoosttag that 

| the asm* procaw wuat opsr* to at Jupisar, 

_ OA* 
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though noc to tht mm ixti nt ThU low of 
He it a likely cKpUnanon of the M« d«pU- 
tkm observed by thi CPMS, if tfc* tolubility 
of N« In He ii m krpe m km been m ifpomH 

(15), 

The vilue of (5 ± 2) X 10 -1 deduced 
for D/H from HQ/Hi • hiriatr than the 
catk* derived horn esrdwbeied ipectro- 
icopfc obMrvAdom (J6. I7h U s 1> x 
10 1 from CHjRfCH* end (2 S l)x tO' 1 
hum HtVHj. ft egrvee well with the value 
(bund in a recent tnalyrii of Voyager infra- 
red ipeetn (J5)i (M * 0.5) X 10" 5 . The 
jovtaa CVH ratio rhouU be the wme nn rhe 
proemnltr or lolnr-nobulrt value because 
the nebula w at tha toura o f Jupiter 1 * by- 
^ofe n . [The tentitrial ratio ia 1.6 X 10" .1 
The tun convened all to 0 to 'He early in 
its binary; thui, Cieta (8) argued that the 
original D/H tn the sola nebula ihould be 
equal to die difference betwmni die 'He/ 
*He ratio* in the tolar wind wwl In meteor- 
ite*, which ihould eon tain nebular He (19). 
Oeiit'i value for pmrrwolar D/H U (2.6 — 
1.0) X 10“*. and thin overlap* the OPMS 
value. Substituting the new determination 
of procosolar ’HePHe ■ (l.l £ 0.2) X 
10“* for tha mecaorltk one. the wine argu- 
ment would yield D/H ■ (3 S l) x 10* *, 
e value that la cioect to that derived directly 
from HD. It h dear that D/H on Jupiter ia 
distinctly hither than the value of (1.6 t 
0.1) X 10“^ found In local Interstellar H 
(30), m axpactad from the tteady destruc- 
tion of D in tha galaxy by ittllar nucleotyn- 
cheats during the hn 4.S X 10* yean. 

Earth -baaed imaging (21) and the re- 
markably until number of cloud particle* 
detected by the nephetomacar axperlment 
(32) indicate that me probe may have da- 
reended Into an uuuauiOy clear pert of the 
Jovian atmosphere. An obvious explanation 
It that the raglon awe one of subsidence, 
like moat riser regions an Earth. Our H,0 
abundance value, then, nay be lew then 
the planetary avenge, although this result ia 
unlikely. In a cloud, the rtletivt humidity 
dm Id be just over 100%, end this U con- 
sistent with the comprahenaivt anetyaii of 
VoyagaT ipeccn by Carbon n eL (18). Tha 
corresponding profile of water vapor he* a 
•cals heigh* of about 3 lua. If. for example, 
die relative humidity it reduced by a down- 
draft frooyfu0% to 90%. the amount of 
wbtidenadi required b 3 ln(90/100) km, ot a 
mere r0-3 Inn. Our measurement! were 
mode tt hr greater depth* than tha expect- 
ed condensation level (3 or 6 bun); to 
perturb the humidity at luch a deep level 
would require • JuwndnJi extending over a 
range of well over a scale height, for more 
then needed limply to clear out the cloud*. 
That* remain* me feasibility of a global- 
acale circulation, perhaps upward at high 
latitudes and downwwd it low. 

The larga atmospheric abundance of C 


suggest* that wurcct other than tha am of 
tht tolar nebula contributed a dgnlneant 
that* of die volatile* In the envelope of 
Jupiter. Those siemens and ewer should 
have been in the eondeneed phew in the 
neighborhood of die accreting Jupiter. Ow 
possible explanation for a deficiency of as- 
ter was proposed by Utuciar and Owen (3, 
4): moat of the plane retires li accreted to 
tha rocky, icy core of the embryo pi* net 
rather than baing collected after die planet 
began to amass its gaseous envelope from 
tha neighboring tolar nebula. Water vapor, 
being me Wait volatile, remained cloae to 
the core, whereas gaits tike CO, CH*, and 
N, could mix with the nebula -dan red en- 
velope, enriching It in C end N. Although 
this scenario could reproduce tha observed 
abundance!, tt mutt be tested by rigorous 
modeling. If there waa a targe mare uf late- 
accreting planctasimels, the lew relative 
abundance of aonosphatie H,0 require* 
them to have been rich in earhonacanue 
materials and deficient in O. It hat been 
customary to iwume diet outer toler-eymm 
plane teiimak reatatble the lew-altered 
carbonaceous chondrites that foil to Earth, 
which contain more O than C (23). but 
there are many that an much drier. Tha 
mean denaiefea observed for many tettUltei 
in the outer tolar system hove been ex- 
plained by awumtag e rock) Ice redo of 
order untry. but perhaps much of the low- 
density material it carbonaceous rather 
than icy. 

If our results are representative of tha 
entire planet (24), they preclude modal* in 
which several Barth manes of late-secret- 
ing watar-nch pknetetlmali deliver vola- 
tiles to the juvlan envelope ( 1 ). There it no 
sign of die dense water eloud, or die corre- 
sponding vapor, invoked to explain waves 
spreading from tire Shoemaker-Levy 9 Im- 
pacts (23). Water vapor observed in the 
Impact plumes probably cane from the lm- 
pectoas, not Jupiter, as lugpated {36). 

MPtMNOB AND NOTH 
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Infrared remote soundings and in situ mass spectrometric measurements place the present water vapor abundance in 
the lower atmosphere of Venus at 30 ± lOppm, and there is some evidence of spatial and temporal variability. 
These measurements also agree that the D/H ratio in this water vapor is 150 ± 30 times that of terrestrial water. No 
other hydrogen compounds are present with mixing ratios as large as 5 ppm v/v. The significance of these 
measurements for the evolution of water in the planet depends on the history of escape of hydrogen and deuterium 
and their replenishment by exogenous and endogenous sources - such as comets and volcanoes. A charge separation 
electric field driven upward flow from the night side bulge region produces escape fluxes averaged over the entire 
planet of (1.5 ± 0.2) x 10 7 cm~ 2 s -1 in H* and (5.6 ± 1) x 10 4 cm -J s -1 in D + during solar maximum (fractionation 
factor - 0.15). Very serious depletion of H* andD* in the bulge during periods of low solar activity reduces the 
escape to 0.15 these rates. The solar cycle average loss rates are (9 ± 1.0) x 10 6 cm‘ 2 s" 1 for H and 
(3.2 ± 0.6) x 10 4 cirr 2 s for D. Until charge exchange loss rates are calculated for atmospheric models that conform 
to Pioneer Venus observations they can only be estimated. The best guess is that H escapes at rates between 9 and 
7 x lO^cnrV* during the solar cycle so that the solar cycle average is 8 x 10*cm -2 s -l .Tbe fractionation factor 
should be small - about 0.02. The total flux could be as large as 1.7 x 10 7 cm~ J s _l and as small as 9 x 10 4 cm“ 2 s _l 
depending on the actual charge exchange contribution. The loss rates can also be calculated from the solar cycle 
modulation of bulge densities. For hydrogen, a solar cycle planet-wide average of (6- 7.5) x 10* cm-^" 1 is 
deduced, f is very large - 0.44. This result suggests that the electric field driven flow in the thermosphere is 
responsible for the entire escape flux of about 8xl0 6 cm’V l ,0.15^f £0.44. If today’s water vapor is 
predominantly a remnant of an ancient supply of water and the fractionation factor has not changed, Venus has lost 
the equivalent of at least 4 meters of liquid water spread uniformly over its surface - 0.12 percent of a full terrestrial 
ocean (fto). This lower limit would be 30 times larger if f is 0.44. If there is a source of water with 
D / H = 1.6 x 10-\ the amount of early water required to account for today's large D/H ratio increases to something 
between 0.6 to 16% of a fto. It is also possible that the water present today has been supplied by outgassing of 
highly fractionated water from the mantle, perhaps during a massive resurfacing event 


1. INTRODUCTION 

Today, the surface and atmosphere of Venus are very dry. How dry is a question that this 
chapter will address. How dry the interior may be is a question much more difficult to answer 
and thus one this chapter will consider only in passing. The planet is losing water because 




hydrogen escapes into space after water molecules have been dissociated in the upper 
atmosphere. Observations made by instruments on the Pioneer Venus Orbiter (PVO) have 
greatly improved our understanding of escape processes and their magnitude. This chapter will 
examine these findings in considerable detail. The planet’s atmosphere may also be acquiring 
water at a rate comparable to the loss rate because of encounters with comets and 
contemporaneous outgassing of the interior. Thus, in principle, an appreciable component of the 
water vapor in the atmosphere today may have been contributed by one or both of these sources. 
On the other hand, much of the water vapor, and particularly the deuterated water vapor, may be 
the residue of water Venus acquired as it was being created. It is also possible that the 
provenance of much of the deuterated water vapor (HDO) is very different than that of ordinary 
vapor (HjO) since the light isotope of hydrogen might be in a virtual steady state balancing 
escape and injection of deuterium-poor water while the heavy isotope is mostly a still decaying 
remnant of the primordial planetary water inventory. In such a case most of the light isotope in 
the atmosphere today could have arrived there relatively recently. This chapter will devote much 
of its attention to a discussion of the degree to which the source of today’s water can be 
determined. In this regard explaining the remarkably high value of the deuterium to hydrogen 
ratio in the water vapor will be the principal objective of the chapter. This review will elaborate 
on developments since the last major reviews on this topic were published (Hunten et ai, 1989, 
Kasting and Toon, 1989). 

2. PLANETOGENESIS: INITIAL VOLATILE INVENTORY 

There is almost a consensus today that the terrestrial planets were formed by collisions 
among planetesimals that grew by accreting gases and dust from the solar nebula during several 
tens of million years after the sun itself developed. Models of the latter stages of this process, 
after 10 23 - 10 26 gm embryos had been formed, place the maturing planetesimals on very 
eccentric orbits that would have thoroughly mixed the material out of which the four planets 
were made (Wetherill, 1991). Some or most of the volatile inventory of the planets, including 
water, may have been acquired during this phase of accretion. In this case, Venus and Earth, in 
particular, should have started with about the same endowment of volatiles, although they may 
have had very divergent experiences with outgassing and release of volatiles from their interiors 
after they reached their present size. Another very important source of volatiles for these planets 
may have been a heavy bombardment by volatile-rich comet-like objects arriving from the outer 
solar system toward the end of the accretionary period. The crater record clearly shows that such 
a bombardment occurred. How rich the intruders were in volatiles is less evident, but it is 
plausible to postulate that they were an important source of volatiles. Again, it is difficult to 
understand how they would have contributed very different amounts of volatiles to Earth and 
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Venus. On the other hand, singular events such as giant impacts by Mars-sized objects that 
caused large scale atmospheric erosion may well have affected one planet and not another. A 
recent attempt to explain the great difference in the rare gas abundances in the atmospheres of 
Earth and Venus (Pepin, 1991) invokes such a giant impact to power an atmospheric blow off 
event on Earth that did not occur on Venus, destroying an initial balance in the two rare gas 
planetary abundances. This model also invokes early blow off of a dense hydrogen atmosphere 
driven by strong extreme ultraviolet radiation from the young sun to explain the elemental and 
isotopic rare gas fractionation patterns. 

An alternative to fractionation by hydrodynamic escape to explain volatile abundances is 
that the terrestrial planets received their argon, krypton and xenon, already fractionated, from 
impacting icy planetesimals (Owen, Bar-Nun and Kleinfeld, 1992). To account for the large 
abundance of argon on Venus, a late impact by a massive comet would be invoked. This 
proposal has been criticized on the ground that the gases would be eroded in the very impact that 
is meant to supply them (Kaula, 1994). A criterion is needed to determine whether irapactors are 
to be sources or eroders of planetary volatiles. 

Further support for the twin planet” model is the near equality of the nitrogen and oxidized 
carbon inventory. This balance is attained by comparing the CO 2 converted to carbonate on 
Earth with atmospheric CCK on Venus. The conversion on Earth depends on water for the 
weathering of silicate rocks and on an ocean in which the carbonate can precipitate. Thus, the 
absence of liquid water on Venus is crucial for the difference in presentation of carbon but begs 
the question of whether there was ever abundant water on that planet Among the volatiles 
whose abundances have been assessed water is the odd one ouL 

On the other hand, Lewis (1972) has shown that during the time when dust and gas were 
condensing to form the primitive solar system condensates, the temperature in the planetary disk 
at 0.7 AU was so high that minerals containing water of hydration would not have been formed 
there. If the terrestrial planets accreted most of their components from material that condensed 
near their present orbits, primitive Earth might have contained much more water than primitive 
Venus did. Indeed, Lewis (1972) argues that the density differences among the terrestrial planets 
can be explained if this is the case, whereas it is difficult to understand the apparently systematic 
density variation otherwise. Prinn and Fegley (1989), however, have pointed out that gas-solid 
reactions are kinetically inhibited at typical nebula temperatures, so that neither Venus nor Earth 
or Mars could have obtained their water this way. 

If, on the other hand. Earth received an appreciable contribution of water from a 
bombardment of volatile-rich objects at the end of accretion, it is hard to understand why Venus 
would not. However, volatiles acquired as a late accreting veneer and those stored in the interior 
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of the planet may well have undergone very different evolutionary paths (Pepin, 1995). In 
particular, the late arriving volatiles may have been severely eroded by planetesimal impact and 
enhanced euv heating not long after their arrival, while those in the interior had to wait for their 
release by outgassing processes. Thus, if for some reason the orbits of large planetesimals did 
not commonly cross each other in the later stages of accretion, the volatile endowment of Earth 
and Venus might have been very different 


3. EVOLUTION OF ESCAPING D AND H 

The measurement long considered diagnostic to the determination of whether or not Venus 
once had abundant water interacting with its atmosphere is that of the ratio of deuterium to 
hydrogen ( D/H ratio) in its water today. Most processes by which hydrogen isotopes escape 
from a planet discriminate strongly against loss of deuterium, leading to enrichment of 
deuterium. If the relative efficiency of deuterium and hydrogen loss -- the fractionation factor -- 
is known, the ratio of the hydrogen content of a reservoir supplying the escaping hydrogen at two 
arbitrary times can be determined from the ratio of the D to H ratios. This is irrespective of the 
absolute magnitudes of the content of the reservoirs or the escape flux (unless the net loss of 
hydrogen has a zero rate). Thus, if the abundance of hydrogen in the atmosphere of Venus today 
(probably in the form of H 2 0 and perhaps some H 2 ) can also be determined, the present D to H 
ratio measured and the D to H ratio in primitive planetary water measured or plausibly estimated, 
the size of the early reservoir can be fixed. 


If the amount of atomic hydrogen in a vertical column of hydrogen compounds is [H] and 
that of deuterium [D] , the respective escape fluxes <pi and 92 , and source strengths P, and P 2 . 


at 

421 

The fractionation factor is normally defined as 

f _ <P2/[P] 

q»/[H] 



( 1 ) 

( 2 ) 

(3) 

(4) 


where 


[ 4 ] 




R = [D]/[H], 

and, in case water is the dominant hydrogen component, 

r= IM22I 

2[H:0 ]' 


(5) 

(5) 


Integration of (4) gives the Rayleigh distillationship relationship for the ratio of the size of the 
reservoir at some time t 2 to a later time r, , 


[H(fi )] '-Rd,/ 


( 6 ) 


independent of <Pi and Pi, if f is constant. If, at r,, [H] is indistinguishable experimentally from 


the hydrogen steady state, [H ss ], that is if <p, = P,, f would be undefined. Eq. 6 then could not 
be used to relate [H] at a very early time, t 2 , to [H] at r, . 

If the effective D to H ratio in the sources is R s , Eq (2) becomes 


d[D] 

dt 


= Rs Pi -Rf (pi. 


(7) 


Eqs. (1) and (7) can be integrated, numerically if necessary, provided there is a way to get a 
handle on the escape rates and source strengths and how they evolve. An analytical solution can 

be found if the hydrogen loss rate is a linear function of the total amount of hydrogen in the 
system, 


<p . = mi (8) 

and both Pj and P 2 are time independent. In that case, 

[H] = [H ss ] + [H 0 - Hss] c'**. 

[D] = [Dss] + [D 0 =-Dss]e- 1K ‘, 

where 

[Hss] = Pi/K = (piss/K, ( 11 ) 

with a similar definition for[D ss ]. [H 0 ] and [Do] are the initial abundances. The ratio of (10) 
and (9) (Gurwell, 1995) gives 


(9) 

( 10 ) 
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( 12 ) 


= Ro([Hq1/Hss) e- ft / T H +(R s /f)(l - e - ft / T " ) 
1 ' ([Ho]/[Hss])e''/™ +(l-e _,/TH ) 


where 


Xh=Hss/9ss. (13) 

In a typical case a hydrogen inventory of 8 x 10 22 atoms per cm 2 is being evacuated at a rate of 
10 7 atoms per cm 2 s so that Th is about 0.25 Gyr. R(r) is a function which increases with time 
to a maximum if Rs is small. It then decreases with time to a limit 


Rss = Rs/f. 


(14) 


This will be the approximate value of R in the extreme or "mature" steady state where neither 
[H] nor [D] is any longer varying appreciably with time. 

[H], of course, reaches an approximate steady state long before [D] does. During the time 
when [H] is virtually in a steady state there is a useful approximation to (12), obtained by setting 
the denominator in (12) equal to unity, or setting 


in (9). Then 


where 


P= K[H] 

(15) 

4 H Lo 

(16) 

dt 


[H] = [Hss] 

(17) 

Ri(t) = Rs/f + [roRo -(Rs/f] 

(18) 

ro = r(o,oo) = [Ho]/[Hss]. 

(19) 


This approximation has appeared several times in the literature with [Hss] set equal to [Ho]. 
Sometimes this has not been appropriate, particularly in the late stages of Rayleigh fractionation 
of an early reservoir [H 0 ], where the hydrogen escape flux and the hydrogen source strengths are 
approximately in balance and Rs very different from R/(r). 

The behavior of R,(t) in the approximation (18) is strongly dependent on r 0 . Long after 
[H] has reached a virtual steady state, [D] can be much larger than [D S s] and R/ very large 
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compare^ to R s /f, its steady state limit. Thus, long after hydrogen has lost all memory of the 
originalldeuterium may retain it. The expression (18) seems to have appeared in the literature 
first in a paper dealing with oxygen fractionation on Mars (McEIroy and Yung, 1976). The 
conditions specified in that paper were that loss of oxygen was balanced with its production from 
water at all times so that from the very beginning oxygen was in the steady state. Enrichment of 
“O occurred because of a fractionation factor effectively determined by diffusive separation of 
l6 0 and u O above the homopause. The counterpart of Eq. (18) was appropriate to that situation 
as it is to the application of Eq. (18) to hydrogen when water is taken always to be in a steady 
state between cometary impacts and escape of hydrogen (Grinspoon, 1987; Grinspoon and 
Lewis, 1988). But it was not to others in which its hydrogen analog with [H 0 ] set equal to [H ss ] 
has been applied to Martian and Cytherean hydrogen. 

In (12) two characteristic times appear. One is the time constant associated with 
establishment of a steady state for hydrogen Xh, the other, for the establishment of a steady state 
for deuterium and, thus, for R is 


XD=T H /f. (20) 

or 1.7 Gyr if f is 0. 15 and x»is 0.25 Gyr. In Figure 1 R(t/T H ) is plotted for several values of r Q 
and f, with R s = R 0 = 1.6 x 10~*. The approximation (18) fits the curves in Fig. 1 for R(i/Xh) 
well if t/tH is greater than 15. In the case of Venus, with a hydrogen escape flux of 10 7 cm~ 2 s _I 
and Hss equal to 8.3 x 10 23 atoms cm’ 2 (corresponding to 30 ppm of H 2 0), r/x would be 15 or 
larger for all times after 4 Gyr. In Figure 1 curves are also plotted for the case 
r oe *^ tH » (1 — e - ^ tH ) corresponding to Rayleigh fractionation of a primitive ocean and ro = 1, 
which is a pure steady state. Practical application of these relationships to the atmosphere of 
Venus will be made in Section 6. 

A useful approximation for r 0 can be obtained by rearranging (12), setting R s = R 0 and 
P = R(t)/R 0 . noting that in all practical cases 

. 5.05 

t/iH> rif (21) 

e _(t/x H ) <<; j (22) 

and 

pe-(^H) « e -f(i/t H ) ^23) 

for 
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tjt H - 8 


(24) 


and any reasonable value of f. Thus, 


r D = pe f,/TH - (e f,/t » - l)/f . 


(25) 


In most cases. 


r 0 =(p-f-‘)e^H, (26) 

since usually 

e ( 'A« » i_ (27) 

r Q is plotted in Fig. 2 along with the approximation (26) for a large range of fractionation factors 
and for p = 150. Such a large fractionation factor cannot be achieved at times 

r<5.05x„/(l-f) (28) 

because Eq. (25) becomes singular at this time. But as t increases above the limit set by Eq. (16) 
the value of r 0 required to achieve p equal to 150 in the presence of a source with 
R s = 1.6 xlO -4 decreases rapidly to a minimum. It then increases exponentially for values of 
r/x H larger than 6 to 8 depending on f. For times much greater than x H , large values of r Q are 
required if an enhancement of D/H as large as 150 is still realized after a long period of quasi- 
steady state between escape and a source of low R s . It is most interesting that, for f > 0.02, an 
enhancement, R, of 150 cannot be obtained at any time unless r 0 is larger than about 180. The 
signature of a massive loss of hydrogen will have persisted, no matter what the hydrogen steady 
state lifetime, if R is large. 

On the other hand, if the planet never had a large amount of water exposed to its atmosphere 
and the sources always have been balanced with escaping hydrogen, so that r D is identically 
unity, R(r) would have risen from R s to R s /f monotonically. An enrichment of 150 times r Q 
would have been attainable with a source R s of 1.6 x 10" 4 only if f is less than 6.4 x 10 -3 . X& 
would then be 50 Gyr. If f is as large as 0.15, R s would have to be 3.75 x 10' 3 , 22 times r 0 , to 
produce a steady state of 150 times R 0 if Ro = 1.6 x 10 -4 . 

Unfortunately, it is not clear that the ideal case in which <pi varies linearly with [H] is often 
attained on any planet including Venus. One practical problem is that the tropospheric hydrogen 
may be mostly bound in condensables such as water vapor, which is prevented by a cold trap 
from passing freely through the tropopause. This is certainly the case on Earth, where 
stratospheric CH 4 and H 2 are the sources of much of the hydrogen above the tropopause. Only 


[ 8 ] 




if the CH 4 and H 2 abundances can be shown to vary linearly with the water vapor content of the 
troposphere will q>i be proportional to [H]. This discussion neglects the water in liquid form, 
which provides a source for atmospheric hydrogen that surely will balance escape until the 
oceans run dry. 

On Venus the clouds constitute a barrier between the lower atmosphere and the upper 
atmosphere. The hydrogen mixing ratio is much less above them than it is below them. It 
remains to be determined how the mixing ratio in the upper atmosphere determining the escape 
flux, would be affected if the amount of water vapor in the lower atmosphere were very much -- 
say orders of magnitude -- larger than it is today. 

In the case of Mars, H 2 0 is optically thick to the solar ultraviolet radiation which 
photolyzes it. Increasing the H 2 0 mixing ratio in the troposphere merely increases the altitude 
at which H is produced from H 2 0. The consequence is that stratospheric H and H 2 , which are 
generated by a complex of chemical reactions involving odd hydrogen and odd oxygen, change 
only slightly, even when H 2 0 changes by very large amounts. Consequently, the hydrogen 
escape flux changes non-linearly with the total amount of hydrogen below the diffusion 
bottleneck near the homopause. 

4. DIAGNOSTICS: REQUIRED 

The properties of the atmosphere to be determined if the problem of the evolution of water is 
to be addressed effectively include 1) the D/H ratio in water or some suitable substance such as 
atomic hydrogen, molecular hydrogen, H 2 S or their ions, 2) the abundance of water and other 
hydrogen constituents of the mixed atmosphere, with their temporal and spatial variations, if any, 
understood and 3) the escape fluxes of hydrogen and deuterium suitably averaged globally, 
diumally, and over the solar cycle. Also required, but more difficult to obtain, are the average 
rates for injfection of H and D by comets and by outgassing of the interior. Thanks mostly to the 
Pioneer Venus mission and remote spectrographic sensing from Earth and the Galileo spacecraft, 
most of the items on the first of these shopping lists by this time have been acquired. It is also 
possible to set reasonable bounds to the rate at which comets inject water, and the D/H ratio for 
water on one comet, P. Halley, is now known with good precision (Balsinger et al., 1995, 
Eberhardt et al., 1995). Magellan has provided enough information about volcanism on Venus to 
allow speculation about the outgassing history of the planet to be at least informed. 

5. DIAGNOSTICS: MEASUREMENTS 
A) D/H Ratio 

By this time the D/H ratio is one of the most thoroughly measured properties of the planet. 
These measurements, all concordant, fall into several classes. 
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1) Measurement of the HD0/H 2 0 ratio in the clouds by mass spectroscopy; 

2) measurement of the HD0/H 2 0 ratio in the mixed atmosphere spectrographically from 
Earth; 

3) measurements of the 1 amu and 2 amu ions in the ionosphere and the interpretation of the 
measurements in terms of ion chemistry; 

4) calculation of the H and D densities in the upper atmosphere from measurements of H + , D + , 
0 + and O densities at altitudes where these ions and neutrals are in a chemical steady state. 

The first of these measurements was made by the neutral mass spectrometer (LNMS) on the 
Large Pioneer Venus Probe as the probe descended. In the clouds, between about 50 km and 25 
km, the inlet leaks to the LNMS were sealed by sulfuric acid droplets. Very large signals in 
several mass channels near 18 amu and 19 amu appeared at this time. The spectrometer was then 
sampling H 2 0 and HDO from the sulfuric acid drops. There were also greatly enhanced signals 
from CH 3 D and HD, apparently created by transfer of deuterium from Cytherean HDO to 
artifact CH* and H 2 in the LNMS (Donahue and Hodges, 1992). Taking account of all the D 
and H represented in HDO, CH 3 D, HD and H 2 0 gave a D/H ratio of (2.5±0.5)x 10* 2 for 
Venus water, 157 ± 30 times the D/H ratio in SMOW. 

McElroy et al. (1982) first pointed out that the ion of mass 2 in the ionosphere could not be 
H 2 + and must, therefore, be D + . They concluded that the measurements implied a hundredfold 
enrichment of deuterium in the atmosphere of Venus, measurement of the densities ions at 2 amu 
and 1 amu by the PVO ion mass spectrometer were analyzed with the help of data from the 
neutral mass spectrometer. The assumption was made on the one hand that the 2 amu species 
was H$ and on the other hand that it was D*. The reactions removing Hj and D* are quite 
different, the first being an ion molecule reaction with C0 2 forming C0 2 H + and H, and the 
second a charge exchange reaction with O forming 0 + and D. The alternatives lead to very 
different height profiles for the 2 amu ions. Clearly the D + option was favored. One set of 
analyses by Hartle and Taylor (1983) gave (2.2±0.6)xl0 -2 for the D/H ratio after 
extrapolation to the homopause. The other set by Kumar and Taylor (1984) gave 1.4 x 10‘ 2 and 
2.5 x 10' 2 for the ratio when corrected for diffusive separation. 

Neutral mass spectrometers have serious problems in measuring the density of neutral 
atomic hydrogen because of its extreme reactivity. Fortunately, the atomic hydrogen and 
deuterium densities in the lower thermosphere may be calculated from measurements of the 
densities of 0 + , H + , D + and O with ion and neutral mass spectrometers in the region where the 
lifetimes of these species for ion molecule reactions are short compared to diffusion times. 
Because the resonant reactions 


H + +0 <-» H + 0 + 


(29) 
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and the H + removal reaction 


H + +C0 2 ->C0H* + 0 (30) 

determine the densities of H + , O and H; in this region the H density may be calculated from the 
relationship 


"(H) = (n(H + )/n(0 + )} {*j9#i(0)(Ti/T)> + *30/»(CO 2 )}, (31) 

where n are number densities, » are the reaction rate coefficient for reactions 29 and 30 and 
T, and T are the ion and neutral temperatures (Brinton, et al., 1980). An analogous expression 
holds for n(D). The measurements of these quantities will be discussed in detail in section (50- 
H and D densities were obtained in the night time bulge region of the thermosphere during the 
first three Cytherean years that PVO was in orbit and during its reentry during 1992. When 
extrapolated to the homopause the ratio of deuterium to hydrogen densities obtained was 
2. 17 x 10~ 2 . Considering the uncertainty introduced by the need to extrapolate the thermospheric 
ratios, the agreement among all these determinations must be considered excellent. 

Near infrared spectral bands of H 2 0 and HDO emitted deep in the atmosphere of Venus in a 
window near 2.3 |im have been analyzed to determine the mixing ratios of H 2 0 and HDO by 
two groups. One working from Mauna Kea (de Bergh et a/,1991) analyzed their observations 
to derive a D/H ratio of (1.9±0.6)x 10' 2 ( 1 20 ± 40 times SMOW) in the mixed atmosphere. 
The other, observing from the Kuiper Observatory Aircraft (Bjoraker, 1992) obtained 
(2.51 0.3) xlO -2 (157118 times SMOW). 

All of these measurements are in good agreement with each other. The ionospheric 
measurement tends to be a little lower than those effectively obtained by observing the mixed 
atmosphere. This may be because of the need to extrapolate measured densities to the 
homopause in the case of the measurements made in the upper atmosphere. Because of the 
possibility of error in this procedure, we shall prefer the “tropospheric” measurements and use a 
value of 2.4 x 10' 2 (150 x SMOW) in this chapter. 

B) Water Vapor Abundance 

There is no shortage of measurements of the water vapor mixing ratio in the atmosphere of 
Venus. Unfortunately, they lack the concordance of the D/H measurement. The measurements 
reported before April, 1992 are recorded in Figure 3. Those obtained by gas chromatographs and 
humidity sensors on Venera 13 and 14 and the Pioneer Venus Large Probe are orders of 
magnitude higher than the rest. Given the difficulties with contamination confronting in situ 
measurements such as these and their gross disagreement with other measurements free of these 
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problems we regard these data as unreliable. Analysis of near infrared spectra of the day sky 
obtained by the Venera 11 and 12 spectrophotometer yields a mixing ratio that decreases 
monotonically from 200 ppm at 50 km to 20 ppm at the surface. This result is perplexing 
because such a profile has been obtained by no other instrument, including other infrared 
soundings on the nightside. Most of the near infrared soundings in the 1.7 and 2.3 Jim H 2 0 
bands radiated on the nightside of Venus and reported before 1992 gave results that clustered 
around 40 ppm from 50 km to low altitudes. There is one exception reported by Bell et al., 
(1991) where one observation gave a mixing ratio of 200 ppm among others much lower. A 
large latitude dependence in water vapor was invoked to account for the variation of infrared net 
fluxes measured by the four Pioneer Venus entry probes (Revercomb et al., 1985). These refer 
to heights above 40 km. The values shown in Fig. 3 decrease with latitude. The lowest are near 
60° , the next near 30° and the very high measurement near the equator. 

The LNMS data at 18 and 19 amu have been analyzed by Donahue and Hodges (1992 and 
1993). Their results are plotted in Fig. 4. The mixing ratio obtained by fitting the H 2 0 counting 
rate profile to the counting rate profile for 36 Ar was 28* 1 ! ppm, between 10 and 26 km, 6.3^ 
ppm at 51.3 km, 4.2^2 PP m at 55.3 km, and an apparent decrease below 10 km to 7± 3 ppm at 
0.9 km. These mixing ratios are normalized to the 36 Ar abundances of 30 ppm obtained by the 
LNMS. The decrease in H 2 0 mixing ratio below 10 km may only be apparent. It is a 
consequence of a failure of the HDO counting rate, which is due mostly to Cytherean water, to 
rise below 10 km as rapidly as the 36 Ar counting rate. The H 2 0 signal is dominated by the 
terrestrial contaminant and could easily mask a variety of profiles for Venus H 2 0. It is quite 
possible that deuterium atoms are being transferred to the very abundant deuterium-poor methane 
in the mass spectrometer. If there is really a decrease in the H 2 0 mixing ratio, there should be 
another hydrogen compound created to take the place of water vapor. There is no evidence for a 
hydrogenic compound in the LNMS spectrum into which the H 2 0 has been transformed by a 
gas phase reaction. At no mass channel belonging to a candidate reaction product can there be 
found evidence for a count rate sufficiently large at 0.9 km to make up the deficiency in 
H 20 there and which then decreases to zero at 10 km. 

Since 1992 the near-infrared spectra of emissions from the dark side of Venus, obtained by 
Crisp et al., (1991) and Bezard et al., (1990), have been simulated by Pollack et al., (1993). 
These data are fit best by a water vapor mixing ratio that is constant at 30 ± lOppm from 40 km 
to about 10 km, but may decrease between 10 km and the surface. The case for the decrease 
below 10 km is by no means robust. The profile obtained by the LNMS (Fig. 4) well represents 
their reconstruction. High-resolution spectra obtained by Bezard et al., (1991) and deBergh et 
al., (1992) in the 1.18 Jim, 1.7 Jim and 2.3 Jim windows are also well fit by a constant mixing 
ratio profile of 30 ±15 ppm. Galileo NIMS measurements also gave 30 ±15 ppm with no 
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horizontal variations exceeding 20% over a wide latitude range (Drossart et al., 1993). Water 
vapor measurements are also discussed by Esposito et al, (this book). 

C) Methane and H 2 

The PVLNMS registered very large signals at 2, 3, 15, 16 and 17 arau (Donahue and 
Hodges, 1993). All of the 2 and 3 amu detected can be accounted for by terrestrial H 2 
emanating from instrument surfaces, probably in the getter pumps and, except when the leaks 
were plugged, from H$ produced in the ion source by electron bombardment of CH 4 . A 
contribution of H 2 from the Venus atmosphere with an appropriately high HD/H 2 ratio 
associated with more than about 2 ppm of H 2 can be excluded, and the data are consistent with 
no H 2 at all. The CH 4 , in addition to that introduced into the LNMS deliberately to help in mass 
peak deconvolution, can be accounted for by methanation, probably occurring on getter surfaces 
and involving instrumental H 2 and CO from the atmosphere. Thus, this CH 4 component 
mimics an atmospheric gas (Donahue and Hodges, 1993). The feature at 17 amu, except when 
the leaks were plugged, and below 10 km, appears to consist of a mixture of 13 C H 4 , and CH5 
generated in the ion source by an ion-molecule reaction involving CHJ and CH 4 . 

D) Hydrogen in the thermosphere: Lyman- a Scattering 

Ultraviolet Spectrometer (OUVS) (Lyman- a spectra from 20 orbits spanning the first three 
Venus years have been analyzed by Paxton et al. (1988). By fitting emission rates to theoretical 
models for resonantly scattered solar Lyman- a they have determined the atomic hydrogen 
number density and vertical flux at the base of the exosphere (200 km). Both remained 
remarkably constant at (6.0±1.5)xl0 4 cm _3 and (7.5±1.5)xl0 7 cm _2 s _1 in the subsolar region 
as solar activity varied. The integrated vertical column of hydrogen above 1 10 km was found to 
be (3.6± 1) x 10 13 cm -2 . Recently, reanalysis of these measurements by Hartle et al. (1995) has 
revised these values to (7.5±2)xl0 4 cm -3 at 200 km and (2.9 ± 0.6) x 10 7 cm _ V. The large 
flux in a region where the escape flux is much smaller than 10 7 cm _2 s"', supplies lateral flow of 
hydrogen to the nightside, where the hydrogen eventually escapes to space or descends below 
100 km and returns to the dayside. Light atoms participate in the planetary-wide circulation of 
winds to the nightside in the thermosphere and back to the dayside in the middle atmosphere. 
Consequently, the densities of the light species H, D and H 2 display a prominent bulge in the 
post-midnight sector of the nighttime hemisphere. The circulation of the high mass major 
atmospheric species is divergence-free. Because of their much greater scale height, this results 
in a net loss of light gases on the dayside in the thermosphere and a net gain on the nightside. 
The consequent nightside bulge in abundance is shifted into the dawn sector by the superrotation 
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of the atmosphere. For further discussion of the effects of superrotation in the upper atmosphere 
is contained in the chapter by Bougher et al (this book). 

E) Escape of Hydrogen and Deuterium 

Until recently, two processes were thought to make important contribution to hydrogen 
escape. Exospheric temperatures are only 300 K on the dayside and 123 K on the nightside. At 
such temperatures, the planet-wide average Jeans escape flux is about 5xl0 5 cm~ 2 s _1 , utterly 
negligible compared to non-thermal loss processes. Collision of energetic oxygen atoms 
produced by dissociative recombination of Oj ions with thermal hydrogen was once considered 
a major source, contributing a planet-wide average flux of 1.2xl0 6 cm -2 s" 1 , with zero 
fractionation factor (Rodriguez et al., 1984) or 3.5 x 10 7 cm" 2 s'* with a large fractionation factor 
of 0.31 according to Gurwell and Yung (1992). However, Hodges (1993) has pointed out that 
these calculations neglected the effect of collisions by the initially hot atoms with thermal 
oxygen that degrade their velocities before they can collide with hydrogen in the exosphere. 
After the effect of these collisions has been accounted for, the hydrogen loss rate following such 
collisions is reduced to insignificance, compared to the dominant processes involving hydrogen 
ions. The first of these is charge exchange in the exosphere between H + ions with greater escape 
velocity in the tail of the hot ion distribution and thermal hydrogen atoms. Two calculation of 
the loss rate due to these collisions gave discordant results. Most of the loss occurs in the 20% of 
the exosphere comprising the bulge on the nightside. The Monte Carlo model of Hodges and 
Tinsley (1986) predicted a planet-wide average loss rate of 2.8xl0 7 cm _2 s _1 for this process. A 
treatment by Rodriguez et al., (1984) gave 3.2 x 10 6 cm“ 2 s -1 . Donahue and Hartle (1992) have 
criticized both models for using unrealistic H and H + distributions. They also pointed out that 
the calculations consider only solar maximum conditions whereas H + ion densities are reduced 
to at most 0.15 their solar maximum values during solar minimum (Figure 5). When solar 
activity is high the ionopause is far from the planet. 0* ions created in the daytime ionosphere 
are transported in great numbers to the nightside, where they create H + ions in the bulge region 
by charge transfer collisions with neutral hydrogen. As solar activity decreases the density of 
ions in the ionosphere decreases and the ionopause closes in. Transport of 0 + to the nightside is 
reduced and production of H + with the bulge consequently decreases considerably. Another 
theoretical effort to recalculate the charge exchange loss rate using the densities measured by 
PVO for a complete range of solar activity is required. In the meanwhile, Donahue and Hartle 
(1992) recommend a rate of 9 x 10 6 cm' 2 s _I for the planet-wide average escape flux during Solar 
Maximum. This recommendation is based on an estimate of the effect of adjusting the properties 
of the two exospheric models in question to realistic ones. Deuterium also is lost by change 
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exchange. The fractionation factor for the process has been estimated to be 0.02 (Krasnopolsky, 
1985). 

Recently, another important escape mechanism for light ions has been identified. It is the 
analog of the terrestrial polar wind in which light ions are accelerated to escape velocity in the 
charge separation electric field that maintains the major ions and electrons of the ionosphere in a 
stable vertical distribution despite their mass difference. Hartle and Grebowsky (1993, 1995) 
have shown that hydrogen and deuterium ions are flowing upward in these fields with more than 
enough velocity to escape. The outward flux of H + in the bulge during solar maximum is 
7.5 x 10 7 cm -2 s _1 , for this mechanism (E) giving a planet-wide average flux of H + during solar 
maximum of 1.5 x 10 7 cm _2 s _l . The D + planet-wide rate was measured to be 5.6 x 10 4 cm* 2 s' 1 . 
The fractionation factor for the process is, thus, 0.15. 

These escape fluxes will vary with solar activity. During the reentry phase of PVO in the 
fall of 1992 the 10.7 cm radio flux was at a level of 125 x KPWcm^H," 1 . In 1978 and 1979 
when the orbiter was exploring the solar-maximum, the 11.7 cm flux was between 200 and 
250 x 10 22 Wcm^Hj -1 (Figure 6). Since the H + ion densities in 1992 were reduced by a factor of 
0.15 and the electric field loss is proportional to the ion densities, the planet-wide average (E) 
flux of H + was reduced to 2.3 x 10 6 cm -2 s _1 . f remained at 0.15. Thus, the solar cycle average 
electric field driven flux <j>j E is 9 x 10 6 cm“ 2 s -1 with a fractionation factor of 0.15. Because the 
charge exchange loss rate (CE) is proportional to the product of the ion and neutral densities it 
should be reduced from 9x 10 fi cm" 2 s _1 to only 7xl0 6 cm -2 s _1 in 1992. If the electric field 
escape flux and the charge exchange flux are additive and independent, the solar maximum 
planet-wide hydrogen escape flux (E + CE) would be 2.4xl0 7 cm -2 s -1 with a fractionation 
factor of 0.10. The low solar activity flux (E + CE) in 1992 would be 8.9xl0 6 cm -2 s _1 with f 
equal to 0.06. The solar cycle average escape flux would be 0.9 and 0.8 x 10 7 cm _2 s _1 for the,(E) 
and (CE) respectively making the solar cycle average flux (EC + CE) 1.7 xl0 7 cm* 2 s _1 and 0.1 
the effective fractionation factor. These measured and estimated fluxes are listed in Table 1. 

F) Transition from Solar Maximum to Solar Minimum: Flux Calculations 

The Pioneer Venus Orbiter reentered the atmosphere in the autumn of 1992. As it did so its 
instruments, including the neutral and ion mass spectrometer, determined ion and neutral 
densities in the nighttime ionsophere. As solar activity declined the hydrogen escape flow out of 
the bulge diminished, but the horizontal inflow from the dayside hemisphere continued 
undiminished. Consequently, the hydrogen and deuterium content of the bulge grew until the 
downward flux had increased sufficiently to balance the inflow. The densities of hydrogen and 
deuterium during the first three years of the mission and in 1992 have been calculated from the 
orbiter data following the technique discussed (Hartle etal., 1996) and in Section 5a. The results 
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are shown in Fig. 7 and 8. The average hydrogen density in the bulge increased by a factor of 
6.7 and that of deuterium by a factor of 4. That this increase is due to the reduced escape flux 
and not to a change in the circulation from the night to daytime hemispheres suggested by the 
smaller deuterium effect and the absence of any modulation for helium (Fig. 8). (In Fig. 9 we 
show the thermospheric D/H ratios during the four years for which data are available. These 
are the basis for the discussion in Section 5a.) 

The number of hydrogen atoms in the thermosphere and exosphere y is governed by the 
simple relationship 

“ = <Ps A/ 2- (tpo +<Pe)Ae (32) 

dt 

where cp s is the upward flux into the thermosphere on the dayside (area A/2), <po is the return 
flux into the mesosphere at the base of the bulge (area Ab, which is about 20% of the planetary 
area A), and cp e is the escape flux from the bulge. It is possible to calculate the escape flux <t> e 
during solar maximum and <t> e from the bulge during solar minimum in the steady state when 
dyldt = 0 in (30). 


2.5<I>s = <Dd + ^t. (33a) 

2.5(}>s =4 >d+<I>i. (33b) 

<I>d and 4 > d are proportional to the density in the bulge during solar maximum and solar 
minimum respectively, therefore 


<D d = 0.15<t> D (34) 

As was noted in 5A, <t>s is known from PVUVS data (Paxton et al., 1987), as revisited by Hartle 
et al. (1996), to be 2.9 x 10 7 cm _2 s _1 . Unfortunately the same is not true of 4>s, but there is no 
apparent reason that atmospheric circulation which is responsible for the flow of hydrogen 
should depend on solar activity. Thus, it is probably safe to assume that <{>$ and <b s are the same. 
In any event diffusion limitation prevents <J> S from exceeding 1.5 
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where eis the appropriate ratio <J>i / In Table 1 the escape fluxes called for by this 
relationship are contained with the measured elective field flux and estimated charge exchange 
flux discussed in Section 5E. 


Table 1 

< 

>i in 10 7 cm“ 2 s _l 

Process 

e 

Eq 33 

Section 5E 

(<Pe) 

E 

0.15 

6.3-5.7 

7.5 

0.9 

cV 

0.73 

7-6.3 

4.5 

0.8 

E + CE 

0.37 

6.6-5.9 

12 

1.7 1 


Eveij^paying due respect to all the uncertainties in the measurements on which these estimates 
are basedyhe electric field escape process alone appears to give reasonable agreement between 
the measured observed value and the one that explains the solar cycle variations in bulge 
densities. It is difficult to reconcile a very large contribution from H + , H charge exchange with 
these results. In any event a very small planetary average flux of about 1.5 x 10 7 cm -2 s -1 during 
solar maximum and 2 x 10 6 cm -2 s -1 during times of low solar activity ot/tiuiuifitfyi is indicated. iS 

The effective H and D escape fluxes (<p) resulting from this analysis of the bulge are also 
tabulated in Table 2. The “synthesis” flux in the table captures both the measured E flux and the 
flux determined from bulge properties in the case of hydrogen, e is assumed to be 0.15. No 
such synthesis is possible for deuterium. 
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Table 2 



(cpi) Hydrogen in 10 7 cm -2 s" 



— 

E 

CE 

E + CE 


§^■1 


0.15 

0.73 

0.37 

0.15 

0.15 


0.9±0.1 

0.8 

1.6 

0.6-0.75 

0.810.2 

f 

0.15 

0.02 

0.10 

0.44 

0.15-0.44 

(tp 2 ) Deuterium in 10 4 cm -2 s _1 


3.2±0.6 

0.5 

3.7 

8-10 

NA 


A similar analysis for deuterium calls for a solar maximum steady state flux, of about 
1-4 — 1.7 x 10 5 cm -2 s _1 planet wide, (8 - 10 x 10 4 cm -2 s - * solar cycle average) with only electric 
field acceleration important. Combined, these fluxes result in a very large fractionation factor 
of 0.44. It is difficult to conceive of a loss mechanism so highly fractionating in the modem 
atmosphere of Venus.The "observed" loss rate for deuterium ions (Hartle and Grebowsky, 1993) 
was only one third as large -- 5.6xl0 4 cm _2 s _1 . This discrepancy could be removed i/*Pyfor 
deuterium were reduced by a factor of about 3, which can be accomplished only by reducing 
D/H to about L0 x 10“ 2 , given the observed density at 220 km during solar maximum Years 1, 2 
and 3. This would be an unacceptably large correction in the HDO/ H 2 O ratio. 

6. EVOLUTION OF WATER 

We shall first consider two cases based on the fluxes discussed in Section 5E. In one the 
measured charge exchange flux will be added to the electric field flux. In the other it will not. In 
the first case, the average escape flux is 1.7 x lO’cnrrV with a fractionation factor of 0.10. In 
the second, it is (911.0)x 10 6 cm -2 s _1 and f is 0.15. If the present hydrogen abundance [H] is 
set by the 30 ppm of water vapor in the lower atmosphere, the characteristic time required to 
establish a steady state in hydrogen is 


Xh =(160-300) Myr. 


(36) 
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Td . to establish a steady state in deuterium and in R, becomes 

X D = (1.6 -2.1) Gyr. (37) 

Finally, the relationship of Rss to R s is 

Rss = (6.7-10) R s . (38) 

Such large values of f do not allow a highly differentiated "mature" steady state if the source is 
not rich in deuterium. 

One end member scenario for evolutions of water on Venus has the hydrogen and deuterium 
today as remnants of an early larger supply of water. In the case of Rayleigh fractionation, the 
factor r in equation (6) for [H 0 ]/[H] is either 262 or 363, depending on choice of f. Originally 
the water on Venus would have been between 3.9 and 5.4 meters deep if spread uniformly over 
the planet’s surface as a liquid. This is between 0.12 and 0.17 percent of a full terrestrial ocean 
(fto). The mixing ratio of this water in the vapor phase would have been between 1.0 and 1.4 
percent if Venus then had the same dense atmosphere it has now. According to Watson et al., 
(1982), vigorous hydrodynamic escape or "blowoff" of hydrogen, driven by solar extreme 
ultraviolet radiation (euv) heating, would have occurred if the mixing ratio of H 2 in the upper 
atmosphere exceeded 2 percent. Because euv from the early sun was almost surely enhanced, 
this coincidence led Donahue et al., (1982) to suggest that Venus may have started with the 
equivalent of a full terrestrial ocean, either in steam or liquid form. Blowoff, for which the 
fractionation factor was close to zero, would have exhausted all but a trace of this water in 280 
Myr, without changing the D/H ratio appreciably. This concept was elaborated by Kasting and 
Pollack (1983) and Kasting, et al. (1984) who calculated the blowoff escape rate that would have 
been reached in a runaway greenhouse induced when Venus developed an atmosphere in which 
an ocean had been converted to steam. They assumed a radiatively equilibrated atmosphere in 
which there were 90 bars of C0 2 and showed that before the sun reached its present luminosity, 
surface temperatures on Venus would become so high that all the water in an ocean would have 
been converted to steam with a surface pressure of 270 bars. Blowoff would then follow. This 
model was superseded by one in which the effect of convection and condensation was taken into 
account. Initial conditions called for a surface pressure of 1 Bar of an N 2 , 0 2 mixture, one bar 
of H 2 0 and a ocean, liquid because the moist convection left the surface temperature below 
100°C. C0 2 was presumed to be converted to carbonate as on Earth. The stratosphere was wet 
in this “moist greenhouse”. Hydrogen was rapidly lost as in the “runaway greenhouse”. A cold 
trap developed only after nearly all of the water was lost, because of the low pressure of the non- 
condensable gases. Thus, a serious defect of the previous model, in which the development of a 
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cold trap in the presence of 90 bars of C0 2 stopped blowoff with 20 bars of water left was 
avoided. In both models, atomic oxygen was presumed to have been lost, swept along in the 
flood of rapidly flowing hydrogen, and also oxidizing the crust. After the blowoff stopped, the 
remaining oxygen could have been reduced to its present insignificant level by reaction with CO. 
This gas should have accompanied outgassing C0 2 in a much larger relative abundance than 
prevails today. After the water had gone, no medium in which carbonate could form would have 
remained, and C0 2 would have accumulated in the atmosphere to its present level. 

As was discussed in Section 2, the simple Rayleigh fractionation relationship (6) cannot be 
used to determine the initial hydrogen inventory if hydrogen is in or near a steady state, that is, if 
Pi = <pi and d[W\ldt = 0. It is not the amount of hydrogen remaining in the atmosphere [H] that 
is decisive in establishing this situation. It is the size of <pi compared to P|. The fact that <pi is 
so small makes the possibility that it is balanced by a source of hydrogen, particularly from 
comets, very great (Grinspoon, 1987; Grinspoon and Lewis, 1988). The Venus D/H ratio is 
much larger than it would be if there were a mature steady state between escape and a source of 
incoming cometary water. The D/H ratio in Comet P. Halley is 3.16 xlO -4 (Balsinger, et al., 
1995; Eberhardt, et al., 1995). So, if f is 0.15, the steady state D/H ratio on Venus when 
cometary water having such a D/H ratio is injected, becomes 2xlO -3 - less than a tenth of its 
present value. As Gurwell has pointed out, if there is a source with R s equal to 1.6 xlO -4 , 
balancing escape of H, a ratio p of 150 today, when t/x H is between 14.5 and 28 can be attained 
only if r Q is between 1300 and 2300 (Figures 1 and 2). Thus, even more early water would be 
needed than if the cometary source were absent 19.5 to 34.5 meters of water, 5 to 6.3 times as 
much as needed in the case of simple Rayleigh fractionation, would be required. On the other 
hand, an outgassing source, such as suggested by Grinspoon (1993) involving a very depleted 
mantle, would be sufficient to balance the current escape flux. The challenge is to demonstrate 
how it can be sufficiently fractionated. 

A similar analysis based on the fluxes deduced in Section 5H 
~(<Pi) = 7.5 x lC^cm^s^.f = 0.44 - calls for the ratio of the primeval reservoir to the present 
one, r in Eq. 6, to be very large — 7700, if there are no sources. This would require 1 15 meters 
of water, 3.5 percent of a terrestrial ocean. In the case of a hydrogen steady state based on Eq. 
(25), with t/t h equal to 12.5, r becomes 3.5 xlO 4 . The required early reservoir would grow to 
520 m of water, 16% of a full terrestrial ocean. Tss would be only 830 Myr, but the Rs required 
for a steady state with R = 2.4xl0 -2 would grow to 1.1 xlO -3 . (Some quantities determined 
here differ slightly from those quoted by Hartle et al. (1990) because it is assumed to be 150 in 
this chapter but 157 in that paper). 

In fact, the influx of cometary water on Venus probably fails to balance the loss of hydrogen 
to space by about an order of magnitude. Shoemaker and colleagues (Shoemaker and Wolfe, 
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1982, Shoemaker et al. 1^0) estimate the fluence of active cometary mass on Earth, and 
consequently also Venus^to 5 x 10 , 7 gm/yr for the past 100 Myr and L 6 x 10 , 7 gm/yr for the past 
Gyr. Shoemaker et al (1994) also argue that the flux of extinct periodic comets may be 
comparable to the flux of active comets. If so and if 30% of the cometary mass is water, the 
effective or average source of coraetary water P x on Venus is presently between 0.7 and 
1.4 x 10 6 cra - 2 s -1 . The average value of Pi over the past billion years would have been a factor 
of 3 larger. If the hydrogen source from comets is that small the estimates of Ho/H from Eq(12), 
based on steady state assumptions, are too large. Those based on simple Rayleigh fractionation 
(Eq 6 ) should be valid unless they are important endogenous sources of highly fractionated water, 
such as thoe Grinspoon (1993), has discussed. A constant ratio K between the escape flux (pi 
and the hydrogen burden of the atmosphere H, now in the neighborhood of 0.3 Gyr 1 - is 
obviously incompatible with an amplification factor Hq/H of approximately 400 given by Eq( 6 ). 
To obtain such a “small” loss of hydrogen, Th must be about 0.75 Gyr. Alternatively, the 
present state of the atmosphere would need to be very close to a steady state 


This would require that 


H ss = 0.9996H 


(39) 


Pi = 0.9996 <pi 


(40) 


It is much easier to accomodate a constant value of ThU> obtain the amplification ratio H 0 /H 
of 7700 called for by the Rayleigh fractionation relationship ( 6 ) if f is as large as 0.44: 

f p.\ _£!i * c ?v t 

1 — l =7700 h (41) 

l <Pj 

This condition would be satisfied if 91 is 7.2 xlO 6 , as given by the bulge study. 
Pi is 1.4 x 10 6 cm - 2 s _1 and the mixing ratio of water vapor is 40ppm instead of 30ppm. 
Alternatively , for example, it can be satisfied with 30ppm of water vapor, 
91 = 5.5xl0 6 cm _ 2 s _1 and Pi = 1.5 x K^cm^s -1 . 

Proposals for rather different scenarios have been made by Grinspoon (1987, 1993) and 
Grinspoon and Lewis (1988). These generally have involved comets or outgassing of the 
interior, as a consequence of volcanism or massive resurfacing events, as sources of hydrogen. 
Given the large fractionation factor that now appears to characterize escape of hydrogen and 
deuterium, it is not possible for planetary escape to approach a mature steady state in which R is 
as large as it is today if cometary sources of water are as poor in deuterium as Comet Halley or 
terrestrial water. The required value of R, is too large. Grinspoon (1993) has proposed as 
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steady state candidates two examples of possible mantle sources. In one case, the outgassing 
water would be fractionated to within an order of magnitude of the present atmospheric value of 
R, the outgassing rates and hydrogen loss rates would be balanced and R close to its steady state 
value. Grinspoon proposed that the mantle of Venus may have attained its enhanced R, because 
of fractionation in the creation of a severely depleted mantle, or because of retention of the large 
D/H signature in water frozen into the mantle after the planet had undergone an early period of 
highly fractionating loss of hydrogen. Alternatively, Grinspoon suggests that a massive 
resurfacing event within the past billion years released water from the interior. The atmosphere 
was then supposed to have experienced rapid loss of hydrogen in which R was enhanced to a 
value even higher than it has today. Subsequent loss of hydrogen would then have caused R to 
decrease to its present value as it declines toward a steady state in both D and H. The suggested 
evolution of R for these examples is shown in Figure 10 (Grinspoon, 1993). Note that the D/H 
range in this figure is based on the measurement of deBergh, et al. (1991), which is centered at 
120 times terrestrial, rather than 150 times terrestrial. 

7. OUTSTANDING PROBLEMS 

The present hydrogen inventory of Venus and the D/H ratio seem to be well understood, 
although there are some observations which make it out to be much wetter than it usually seems 
to be. Recalculation of charge exchange loss rates, which are based on our present understanding 
of the temporal history of Venus neutral and ionized hydrogen densities and their variation with 
solar activity, are needed. The type of calculation once performed by Kumar, et al. (1983) to 
determine the evolution of hydrogen and deuterium loss mechanisms over the 4.5 Gyr history of 
the planet but based on our present understanding of loss mechanisms, such as the electric field 
process, needs to be repeated. In view of the long time required for deuterium to reach a steady 
state, it would be instructive to investigate the effect on the present D/H ratio and water 
abundance of including the contribution of the remnants of a primitive "ocean" 10 meters deep 
on "steady state" degassing models, such as those proposed by Grinspoon (1993). Models of 
catastrophic mantle outgassing as sources of today’s highly fractionated water vapor need to be 
bolstered with more robust and detailed foundations than those with which they are currendy 
provided. It would be useful to have a study of the Venus thermosphere and mesosphere that 
continuously spans a complete solar cycle. 
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Abstract A strong solar cycle variation in hydrogen and deuterium densities is observed in the 
nightside thermosphere of Venus when Pioneer Venus Orbiter (PVO) measurements made during 
the first three Venus years of the mission are compared with those made during the preentry phase 
of the mission. Solar maximum conditions prevailed during the former period, while near solar 
minimum conditions occurred during the latter. Pronounced density enhancements in H of 6.5 
times and D of 4 times are observed in the nightside bulge region as solar activity decreased from 
maximum to near-minimum values. We attribute the buildup of H and D to a reduction in the 
escape fluxes as solar activity decreased, a behavior that is consistent with the known properties of 
the dominant escape processes due to the charge separation electric field (E) and charge exchange 
(CE). Application of this simple concept leads to expressions for the H and D escape fluxes 
which relate solar cycle variations of the respective bulge densities and escape fluxes to the dayside 
source fluxes. Planet averaged escape fluxes in the ranges 1.23 x 10 7 cm' 2 s' 1 £ 0> tf (H) £ 1.45 x 
10 7 cm' 2 s* 1 and 1.61 x 10 5 cm* 2 s 1 < o, (D) < 2.15 x 10 5 cm' 2 s 1 are obtained by this method 
without specifying any particular escape mechanism. Considering the uncertainties in the 
measured parameters, these flux ranges are in reasonable agreement with the escape fluxes derived 
from the specific processes, E and CE. When fractionation without or with an external source of 
water (e.g., comets) is applied to the range of possible escape fluxes, a fractionation factor /= 

0.44 is obtained and an ancient reservoir of water equivalent to the range 125 m to 570 m of liquid 
uniformly distributed on the surface or 3.7% to 17% of a full terrestrial ocean is derived. When 
the specific escape processes E and E + CE are considered, the /values are 0.15 and 0.1 and the 
magnitudes of the reservoirs are lower, having a range of equivalent depths of 4.1 m to 36.8 m, 
respectively. 


1. Introduction 

Knowledge of the density distribution of hydrogen and its 
isotope deuterium in the atmosphere of Venus is essential to 
understanding the evolution of these constituents and that of 
water. Measurement of H and D densities and determination of 
their escape rates as well as the D/H ratio were some of the 
primary objectives of the Pioneer Venus (PV) mission which 
began exploration of Venus in December 1978. The 
objectives of this paper are to determine the diurnal variation 
of H and D at high and low solar activity, evaluate the solar 
cycle variations of the respective escape fluxes, and apply 
them to an estimate of the size of an earlier water reservoir. A 
new approach is used here to evaluate the escape fluxes which 
requires little or no knowledge of the actual escape 
mechanisms. 
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Early measurements of H were obtained from UV 
instruments on Mariner 5 [Barth et a/., 1967] and Veneras 10 
and 11 [Bertaux et ai, 1982], while more detailed knowledge 
of its density and structure was obtained by several instruments 
on board the Pioneer Venus Orbiter (PVO). Using ion densities 
of IT and 0 + measured by the Orbiter ion mass spectrometer 
(OIMS) and neutral densities of O and CO, observed by the 
Orbiter neutral mass spectrometer (ONMS), Brinton et ai 
[1980] obtained H densities in the chemical equilibrium region 
spanning the first diurnal cycle of the PVO mission. They 
discovered a large diurnal variation of H having a nighttime 
density bulge peaking near 0300 UT which was some 400 
times more dense than that observed on the dayside. Taylor et 
ai [1985] extended the analysis of the H bulge region to 
include the first three Venus years of the PVO. An analysis by 
Paxton et ai [1988] of the Lyman a fluxes measured by the 
Orbiter ultraviolet spectrometer (OUVS) produced H 
concentrations similar to those of Brinton et ai [1980]. 

The first detection and measurement of D in the atmosphere 
of Venus was confirmed after considerable analysis of two sets 


4525 



4526 


HARTLEETAL: SOLAR CYCLE VARIATIONS AND ESCAPE ON VENUS 


of measurements made during the PVO mission. For one set, a 
detailed study of mass spectra obtained from hydrated liquid 
cloud drops that fortuitously clogged the inlet leak system on 
the Pioneer Venus large probe neutral mass spectrometer 
(LNMS) led to the identification of D and an estimate of the £ 

D/H ratio of (1.6 ± 0.2) x 10 2 in the lower atmosphere * 

[Donahue et al , 1982]. At about the same time, McElroy et al J 
[1982] suggested that the mass 2 ion observed by the OIMS 
could be D* instead of H 2 + . Subsequently, an analysis of the 
height distribution of the mass 2 ion measured in the chemical 
equilibrium region of the ionosphere by the OIMS led to the 
clear identification of D + , from which a D/H ratio of (2.2 ± 0.6) 
x 10' 2 was derived at the turbopause [Hartle and Taylor , 1983]. 

A recent analysis of the LNMS measurements by Donahue and 
Hodges [1992] revised the D/H ratio upward to (2.5 ± 0.6) x 
10 ' 2 . 

In a recent study, Grebowsky et al. [1995] expanded the 
database used by Brinton et al. [1980] and Taylor et al. [1985] 
to derive the H distributions. The number of data points was 
increased by testing each derived H density to determine 
whether or not chemical equilibrium prevailed at each 
measurement point. In the previous studies, maximum 
altitudes of the chemical equilibrium region were preset, 
thereby eliminating valid H measurements when chemical 
equilibrium existed above such altitudes. The work of 
Grebowsky et al. [1995] included H densities from the first 
three Venus years as well as those from the preentry phase of 
the mission. The current study will extend this work by adding 
D to the H density database, which will then be used to evaluate 
their respective escape rates. In turn, these results are applied 
to an analysis of the evolution of H and D in the atmosphere of 
Venus to provide an estimate of the size of an early water 
reservoir on the planet. 

2 . Solar Activity and Escape Mechanisms 

Periapsis of PVO was in the thermosphere during the first 
three Venus years of the mission before it ascended into the 
exosphere. During this period, between December 1978 and 


Solar Activity 



Figure 1. The 10.7 cm solar radio flux line spanning the 
calendar years of the PVO mission. The time intervals depicted 
by the vertical lines labled 1, 2, 3, and P correspond to the 
periods when periapsis of PVO was in the midnight to dawn 
sector (only to 0430 UT for P) during Venus years 1, 2, and 3 , 
and the preentry period, respectivly. 



Figure 2. Height profiles of H\ D + , and O* for solar 
maximum (dashed lines) and near-solar minimum (solid lines). 


July 1980, solar activity was at maximum levels, as indicated 
by the 10.7-cm solar flux line shown in Figure 1. In the 
preentry phase of the mission in 1992, almost 14 years after 
initial entry into the atmosphere, periapsis of PVO returned 
briefly to the thermosphere, where measurements were made in 
the predawn H bulge region, just before bumup of the orbiter 
on October 8, 1992. Solar activity was near a minimum during 
this phase of the mission, making it possible to study solar 
cycle effects on escape and evolution when this data set is 
considered together with that of 1978-1980. 

Immediately following the low solar activity measurements 
of 1992, the first study of this kind was done by Donahue and 
Hartle [1992], who noted a strong solar cycle variation in ion 
composition requiring changes in the known H and D escape 
fluxes. An example of ion composition dependence on solar 
activity is shown in Figure 2, where vertical profiles of 
average H\ D*, and 0 + densities in the 0200-0400 UT dawn 
sector are plotted. Two sets of ion densities are shown in the 
figure, one corresponding to averages taken over the first three 
Venus years and the other over the preentry phase, when solar 
maximum and near solar minimum conditions prevailed, 
respectively. The dawn sector is shown here because this is 
the bulge region where H* and D*, and H and D maximize, where 
most of these atoms and ions in the thermosphere reside, and 
consequently, where their respective escape fluxes maximize. 
In this sector, H* and D* densities near solar minimum are 
about 7 times and 3 times less, respectively, than those at 
solar maximum while the corresponding CT densities are a 
factor of 20 times less at low altitudes (160-200 km) with the 
factor decreasing to 7 or 8 at higher altitudes. Similar solar 
cycle variations were described by Kar et al. [1994]. The 
reduction in ionospheric densities on the nightside is expected 
as solar activity declines. This is due to at least two major 
factors. (1) Fewer ions are produced on the dayside as solar 
activity decreases; therefore fewer ions are available to be 
transported to the nightside. (2) With a reduced ion population 
on the dayside, the ionopause altitude decreases, thereby 
reducing the area of the terminator annulus through which day 
to night ion flow takes place. 

The discovery of the strong decrease in FT and D + densities 
as solar activity decreased led Donahue and Hartle [1992] to 
reevaluate the known escape mechanisms. Four basic 
mechanisms were known at the time, two of which are not 
considered to be important today. Thermal or Jeans escape was 
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shown to be insignificant due to the very low exospheric 
temperatures of 284°K at noon and 123°K in the nightside H 
bulge of Venus (note: the temperatures used throughout are 
from the Venus international reference atmosphere (VIRA) 
model [Keating et a/., 1985]). Collisions of H and D with 
energetic oxygen, O*, produced by dissociative recombination 
of 0 2 *, was considered to be an important escape mechanism 
for a number of years [Rodriguiz et al, 1984; Gurwell and 
Young , 1993]. However, Hodges [1993] has recently shown 
that when proper account is taken of the scattering down in 
energy of hot oxygen colliding with ambient thermal oxygen, 
the energy density of the hot oxygen population is reduced to 
the point that the resulting H and D escape fluxes are rendered 
unimportant. This leaves two escape mechanisms that are 
thought to remain important. One is due to charge exchange 
(CE) of “hot” ionospheric H* with H, which produces a 
population of “hot” H* atoms possessing a significant number 
of escape trajectories. The most recent calculations of the loss 
rates due to CE were made by Rodriguiz et al. [1984] and 
Hodges and Tinsley [1982], who obtained planet average 
escape rates of 3.2 x 10* cm' 2 s' 1 and of 2.8 x 10 7 cm' 2 s \ 
respectively. Donahue and Hartle [1992] have shown that the 
large differences in these escape fluxes are due to the fact that 
both investigator groups used widely differing H* and H 
densities, none of which were representative of average 
conditions. When average IT and H densities derived from 
PVO data were used, Donahue and Hartle [1992] obtained an H 
escape flux due to CE of 0.9 x 10 7 cm 2 s' 1 for the solar 
maximum conditions of the first three Venus years of the PVO 
mission. The CE escape process for D is very inefficient 
relative to that for H, and, consequently its fractionation factor 
is only 0.02 [Krasnopolsky, 1985]. The remaining escape 
mechanism of significance is a fluid-like upward flow of H* and 
D\ which is accelerated to escape speeds by the charge 
separation electric field in the ionosphere [Hartle and 
Grebowsky, 1995]. As the plasma flows to higher altitudes in 
the ion exosphere, the motional electric field of the solar wind 
impressed across the draped magnetic field lines of the ionotail 
of Venus eventually overtakes the polarization electric field 
and accelerates the ions up to solar wind speeds as the ionotail 
merges into the interplanetary medium. Essentially all escape 
of H* and D* by the electric field process (E) occurs in the light 
ion bulge, where most of these ions reside. When the escape 
fluxes derived by Hartle and Grebowsky [1993, 1995] are 
averaged over the planet, we obtain values of 1.5 x 10 7 cm' 1 s' 1 
and 5.6 x 10 4 cm’ 2 s' 1 for H* and D*, respectively, for solar 
maximum conditions. 


3. H and D Density Distributions 


Apart from minuscule differences due to their nuclear 
masses, the atomic properties of H and D are the same. 
Consequently, H and D atoms and ions obey the same principal 
chemical reactions 


H + +0 <=> 0 + + H 
D^+0<=>0 + +D 

h + +co 2 «>coh + +o 
d + +co 2 <=>cod + +o 


which determine the H and D densities 


n( H) = 


n(H+) 
n( 0 + J 




(T V /2 


- K 2 n(C 0 2 ) 


(la) 

(lb) 

(lc) 

(l d) 


(2a) 


n( D) = 


n(D+) 

n<0+) 


AT, n( O, 




1/2 


+ K 2 n(C0 2 ) 


(2b) 


at all altitudes in the ionosphere of Venus where chemical 
equilibrium holds. The ion densities n(H*), n(D*), and n(CO 
were measured by the OIMS, the neutral densities /i(0) and 
n(C0 2 ) and their temperature T were measured by the ONMS, 
and the ion temperature 7] was measured by the Orbiter 
retarding potential analyzer (ORPA). The coefficient K l for 
reaction (la) is that used by Brinton et al. [1980]. For reaction 
(lb) we take K{ = K v This can be seen to be a valid 
approximation upon first considering that K x (T t IT )' n can be 
expressed as the ratio of the forward rate factor c f (SkT, /Tip, 16 ) l/2 
to the reverse rate factor a f (8^T/7tp, 16 ) 1/2 near thermal 
equilibrium (as is the case here), where the cross sections ay and 
o f are energy averaged and |i U6 = + m 0 ) is the reduced 

mass of the reactants [e.g., Holier and Banks , 1969], In this 
form, K x = a f la r is independent of mass. Therefore, since K x 
has the same form as K x and the ratio q/a r is essentially the 
same for hydrogen and deuterium collisions, then K x = K v The 
coefficient K 2 is the ratio of the rate 3.8 x 10 9 cm 3 s ' 
[Anicich, 1993] of the reaction H + + C0 2 divided by the rate of 
the reaction O* + H, which is 2.5 x 10'" T m cm 3 s 1 [Nagy et 
al., 1980]. Similarly, the coefficient is the ratio of the 
rate of the D* + C0 2 reaction, 3.5 x 10 9 cm 3 s 1 [Anicich, 
1993], divided by the rate of the O* + D reaction, which has not 
been measured yet. However, considering the simple 
expressions above for the rate factors, the later rate can be 
estimated by multiplying the rate for 0 + + H by (p, J p M6 ) l/2 = 
0.73. 

The principal database used in the analysis is the Pioneer 
Venus unified abstract data system (UADS) consisting of 12-s 
averages of the measurements. When measurements of some of 
the parameters are missing in UADS, empirical models are used 
to fill in the gaps; i.e., the Hedin et al. [1983] model is used 
for the O and CO : densities and the neutral temperature, while 
the VIRA model is used for the ion temperature. In contrast, 
models are not used for the ion densities because they cannot 
properly account for the high spatial and temporal variations 
required by (2). Consequently, when ion densities are missing 
from the database, no H or D densities are derived. 

The chemical equilibrium region where (2a) and (2b) are 
valid is determined by an algorithm developed previously to 
derive H densities [Grebowsky et ai, 1995]. In this case, 
values of n(H) and n(D) are only retained when the chemical 



Local Time (hr) 


Figure 3. Altitude versus local time of each point where 
both chemical equilibrium and a value of n(H) exist. 
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timescale l/*,n(0) is less than the H* diffusion timescale h 2 id , 
where h is the O scale height and d is the diffusion coefficient 
for FT ions [Banks and Kockarts, 1973]. This approach is 
conservative for evaluating n(D) because the diffusion time 
constant for D* is greater than that for H*, and consequently, 
chemical equilibrium for D is valid to somewhat higher 
altitudes than it is for H. Collisions with all major ion and 
neutral species are accounted for in the calculation of the 
diffusion coefficient. The scale height for C0 2 is not used in 
the diffusion time because it is only important at low altitudes, 
where H* is not detected. Applying these guidelines, we find 
that the altitudes where both chemical equilibrium exists and 
enough data exist to allow n(H) to be derived are distributed in 
a distinct diurnal pattern as indicated in Figure 3. In particular, 
chemical equilibrium extends up to about 250-260 km on the 
dayside and only to 170-175 km on the nightside, a diurnal 
variation which is largely a reflection of the variation of the 
major ion CT at these altitudes. The lower limits of the altitude 
distribution in Figure 3 are not the lower limits of the chemical 
equilibrium region, but simply the lowest altitudes where H* 
remained detectable. 

The H and D densities derived from (2a) and (2b) in the 
chemical equilibrium region are shown in Figure 4, where the 
large H density bulge discovered by Brinton et ai [1980] is 
displayed. Helium densities [e g., Kasprzak et ai, 1993; 
Grebowsky et al.> 1995] at 170 km are also shown for 
comparison in Figure 5. An absence of telemetry caused the 
postdusk data gap in Venus year 2, while a special OIMS 
operational mode resulted in a similar gap in Venus year 3 . 
After about 1000 UT in Venus year 3, periapsis of PVO rose 
above the chemical equilibrium region and did not return until 
the preentry phase, where measurements were fortunately made 
in the bulge, just before bumup. Since the latitude of periapsis 
of PVO was about 16° during the first three Venus years and 
about -10° at preentry, the diurnal variations shown in Figure 4 


are considered essentially equatorial for the purposes of this 
paper. The variance in the densities at a fixed local time is due 
to atmospheric waves and structures; some structures have 
been shown by Taylor et ai [1985] to be related to solar wind 
ram pressure changes. Average values of the H and D densities 
are shown in Figure 6, which are convenient for making 
comparisons and carrying out quantitative analysis. 

Apart from the large nightside H and D density bulges, 
perhaps the most striking features in Figure 6 are the 
enhancements of the preentry bulges over those observed 
during the first three Venus years. In particular, we find that 
the density of the preentry H bulge is about 6.5 times greater 
than the average density found in the bulge during the first 
three Venus years. A similar comparison made with D reveals a 
density enhancement of 4 times. Referring to Figure l, we see 
that solar activity levels during the first three Venus years, 
1978-1980, are similar to the solar maximum levels just 
preceding preentry. In the ensuing analysis, we make the 
reasonable assumption that the H and D densities in 1991 were 
similar to those in the 1978-1980 period. In this scenario, 
density enhancements developed when solar activity dropped 
from its high level in early 1992 to the low level of the last 4 
months of the mission. It is fortunate that solar activity was 
relatively steady during this latter period, when all the 
preentry bulge densities were measured, as will be clear in the 
steady state analysis that follows. 

The major gas densities and temperatures in the H and D bulge 
region were observed to be similar during both solar maximum 
and near solar minimum conditions [Kasprzak et ai, 1993]. 
Keating and Hsu [ 1993] have observed that the thermospheric 
temperatures and densities on the dayside decrease somewhat as 
solar activity decreases from maximum to mimimum 
conditions. These observations are consistent with the solar 
cycle variations predicted by Bougher et al. [1986], whose 
circulation model also suggests that thermospheric winds vary 
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Figure 4. Diurnal variations of H and D densities derived from (2) for each Venus year and the preentry period 
at altitudes shown in Figure 3. 
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Figure 5. Diurnal variations of He densities for each 
ONMS in the altitude range 170-180 km. 


only slightly over a solar cycle. The weak solar cycle 
variation in thermospheric circulation, which transports H and 
D into the nightside, can not account for the large solar cycle 
variations of these constituents . However, the weak solar 
cycle variation in thermospheric circulation is consistent with 
the He density observations shown in Figure 5, indicating no 
noticeable solar cycle variation, because the density buildup of 
this species on the nightside is primarily driven by 
thermospheric circulation and not at all by escape [Mengel et 
ai , 1989]. 

4. Escape Rates From Solar Cycle Variations 

We suggest that a likely cause of the H and D enhancements 
occurring in the bulge region as solar activity decreased is a 



Local Time (hrs) 

Figure 6. Diurnal variations of average values of H and D 
densities derived from data in Figure 4 for Venus year 1 (light 
solid line), year 2 (long-dashed line), year 3 (short-dashed 
line), and preentry period (heavy solid line). 
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Venus year and the pre-entry period as measured by 


consequence of the known decrease in the escape rates of these 
constituents [ Donahue and Hartle, 1992] as solar activity 
declines. Thus we examine this concept here and show how 
such solar cycle variations can be used to estimate the escape 
fluxes of H and D. We note that this scenario for H and D is 
fundamentally different from and consistent with that for He, 
which does not exhibit a solar cycle variation and whose 
density variation is not affected by escape, but primarily by 
thermospheric circulation. 

Considering the variation in H first, which can be described 
in terms of the rate of change of the number of H atoms in the 
thermosphere, y, by 

dV 

— = '9s A d~< , ?e + ^d) A b (3) 

at 

where <p 5 is the source flux of atomic H which flows up into the 
thermosphere from the mesosphere after being formed there by 
photolysis of H 2 Q and H 2 . Since the source is confined to the 
dayside, its flux is multiplied by A d , the area of the dayside. 
As H flows from day to night, a fraction escapes from the top 
of the nighttime atmosphere with a flux (p g , while the 
remainder recirculates in the mixing region after flowing 
downward through the lower thermosphere with a flux 
Essentially all of the hydrogen in the thermosphere resides in 
the nighttime bulge region; therefore the downward flow and 
escape flux are multiplied by the area of the bulge, A b , which is 
about 20% of the total area A. We have varying degrees of 
knowledge of y, <p 5 , and cp^ which can be used to estimate the 
escape flux tp^ as will be shown below. 

4.1 Steady State Analysis 

Since solar maximum conditions prevailed just before the 
final entry of PVO’s periapsis into the thermosphere and since 
the corresponding levels of the 10.7-cm solar flux were similar 
to those of the first three Venus years of the mission, the 
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Table 1. Escape Fluxes (10 7 cm' 2 s' 1 ) 


Process 



o. 

o. 

< (p f > 

<?<> 

E 

7.5 

1.5 

1.1 

0.23 

4.3 

0.86 

CE 

4.5 

0.9 

3.3 

0.64 

3.9 

0.77 

E+CE 

12 

2.4 

4.4 

0.85 

8.2 

1.60 


conditions in the thermosphere during the latter period are 
assumed to approximate those of the former period. 
Consequently, when solar maximum occurred during the 

preentry phase, dy/dt = 0, and (3) becomes 

0 s A d =(4> d +Q e )A b (4) 

where <I> refers to solar maximum conditions. Time 

independence also applies to the period when the preentry 
measurements of H were made in 1992 because, as noted 
previously, the solar flux remained steady during that interval. 
Therefore 

= ($d (5) 

where refers to conditions near solar minimum. The steady 
state approximation should be valid over most of the solar 
cycle and not just during the maximum and minimum periods 
considered here. For example, the diffusion timescale h : fk, in 
terms of the scale height h for CO ; and the eddy diffusion 
coefficient k, is only about 0.7 days at 1 10 km and 4 days at 
100 km, periods much shorter than that of a solar cycle. 
Furthermore, the residence time, for thermospheric H is 

only about 10 to 20 days. 

A parameter that will be useful in the analysis is the escape 
flux ratio 

e = <t>, /O,. (6) 

As will be clear below, the relation between and <& d is 
known, <P S is known as are upper and lower bounds on <J>^. 
When these relationships and equations (4) and (5) are 
combined, an expression relating to results. However, 
the magnitudes of <t> e and can be evaluated using current 
knowledge of the dominant escape fluxes [ Donahue and Hartle, 
1992] dbe to charge exchange (CE) and the charge separation 
electric field mechanism (E) to constrain e. Values of the 


Table 2, Escape Efficiencies and Globally Averaged Escape 
Fluxes (10 7 cm' 2 s* 1 ) 


Process 

e 

0 

{From (9)) 

<*>e 

(Table 1) 


0.00 

0.85 

1.23 


E 

0.15 

0.87 

1.26 

1.50 

E + CE 

0.37 

0.90 

131 

2.40 

CE 

0.73 

0.95 

1.38 

090 


1.00 

1.00 

1.45 



escape fluxes obtained by Donahue and Hartle [1992] are 
shown in Table 1. The overbar used throughout corresponds to 
a global average (e g., <t> e - <D^ b M = /5), and the 

quantity <<1>^> = (<P e + <t> tf )/2 is the average flux over a solar 
cycle. The low solar activity escape fluxes were taken to be 
negligible by Donahue and Hartle [1992], whereas they have 
been scaled here by multiplying the solar maximum fluxes by 
n!N* for E and nn!(N*N) for CE. When flux ratios are derived 
from Table 1, e values of 0.15, 0.73, and 0.37 are obtained for 
the processes E, CE, and E + CE, respectively. 

The most recent work pertaining to the source flux has been 
done by Paxton et al. [1988], who derived the upward flux of H 
from Lyman a intensities observed by the OUVS on PVO. 
They obtained an upward flux of 7.5 x 10 7 cm' 2 s' 1 in the 
subsolar region using data from the first three years of the PVO 
mission, when solar maximum conditions prevailed. We have 
reevaluated this flux in the appendix and obtain an upward flux 
<l> s of 2.3 x 10 7 cm 2 s \ which does not exceed the limiting 
flux of <t>, =3.5 x 10 7 cm 2 s' 1 , as required. Furthermore, the 
observed H densities described in section 3 are 1.25 times 
greater than the values used by Paxton et al. [1988]; therefore 
adjusting the upward flux by this factor leads to d> 5 = 2.9 x 10 7 
cm 2 s' 1 . At this time we do not know the value of the source 
flux <t> 5 ; however, it is probably larger than as more H is 
circulated through the upper atmosphere. An upper bound is 
fixed by the fact that the source rate cannot exceed the 
maximum rate Thus there is an upper bound for 0^ = 1,5 

x <t> s . The bounds on the source flux near solar minimum are 
expressed by the flux ratio 

1 ^ " 1-5. (7) 

A downward flux of H in the nightside bulge region has been 
shown to be significant in a number of global circulation 
models [Mayr et el., 1978; Hartle et al., 1978; Mayr et el., 
1980]. To a first approximation, it is proportional to the 
density in the bulge region as expected from the diffusion 
equation. Thus, we define the flux ratio 

*n= <t> d / t>d = N/n {8) 

in terms of the average H densities N and n observed in the H 
bulge during solar maximum and near-solar minimum, 
respectively. When the average densities in Figure 6 are 
scaled, e rt =0.21, 0.17, and 0.09 for Venus years l, 2 and 3. 
respectively, and an average e n = 0,15. These ratios reflect 
density variations in the height range of about 150 km to 170 
km. the limited region where H and D densities could be 
evaluated on the nightside. Consequently, the analysis here is 
restricted to a lower boundary of about 160 km through which 
the source flux and the downward flux flow. The upper 
boundary through which the escape flux flows is taken to be at 
an altitude high enough to include all the significant processes 
contributing to E and CE while avoiding the complicating 
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Table 3. Escape Fluxes (10 7 cm 2 s' 1 ) From (9) 


effects of lateral transport; in particular, we choose an upper 
boundary of 500 km, the approximate altitude of the ion 
exobase which is well above the neutral exobase. 

When values for £, £„, and are used in (4) and (5), the 
escape flux and its time average can be evaluated. The analysis 
is simplified by combining these equations to obtain the 
following convenient forms: 

<t> € =Tife,e n ,8 5 )0 5 (9) 


is the efficiency for escape. Considering the estimated values 
of and e 5 , physical limits can be placed on r\ as e is varied. 
An upper bound J] < 1 results from the physical limit <J> e < 
Q> s , which also requires e < 1 . In the other extreme, when e = 
0, z n = 0.15 (average value) r) = 0.85 and O, = 0.85 . 

Altogether, the following range of variables apply: 0 < e < 1 ; 
0 < q < 1; 0.85 <I >5 < O, < ^ 5 or 1.23 x 10 7 cm ' 2 s ' 1 < < 

1.45 x 10 7 cm ' 2 s' 1 . These limits are applicable when = 1 ; 
however, considering the upper limit z s = 1 . 5 , the ratio r|(e, 
£ n , 1.5)/ r|(e, e n , 1.0) = 0.91 indicates that an ignorable 
change in the fluxes occurs as e 5 spans its limits.. Escape 
efficiencies and globally averaged escape fluxes from ( 9 ) are 
shown in Table 2 for the limiting values of £ and those of the 
three processes E, CE, and E + CE, when e 5 = 1 . Also shown i n 
Table 2, for comparison, are escape fluxes from Table 1 . 

Shown in Table 3 are all types of escape fluxes obtained 
from (9) for the processes E, CE, and E + CE, for completeness. 

A few things are immediately apparent from the range of 
parameters discussed above and in Tables 2 and 3. (1) The 

escape flux in Table 1 which combines the electric field and 
charge exchange processes <t> f (E + CE) = 2.4 x 10 7 cm 2 s 1 is 
considerably larger than the upper limit for the escape fluxes 
inferred from the solar cycle variations in the bulge. (2) 
Separately, the escape fluxes due to the charge separation 
electric field, 4>,(E), and charge exchange, d^fCE), are also 


outside, but very close to the upper and lower bounds, 
respectively, allowed by (9). (3) The escape efficiencies and 
the escape fluxes are constrained to vary a maximum of only 
18% when 0 ^ £ < 1, and no more than 9% when the processes 
E, CE, and E + CE are considered. It is fortunate that the escape 
flux is so tightly bound to the range 1.23 x 10 7 cm 2 s 1 < <!>,< 
1.45 x I0 7 cm' 2 s'\ making it possible to study the evolution 
of H in this flux range without identifying any particular 
escape mechanism. 

Equations (6) through (10) can be applied directly to 
evaluate escape parameters for D. Density ratios z n (D) = 0.25, 
0.4, and 0.1 are obtained by scaling the average D densities in 
Figure 6 for Venus years 1, 2, and 3, respectively, yielding an 
average £ n (D) = 0.25. The physically derived escape flux due 
to electric field acceleration at solar maximum [Hartle and 
Grebowsky, 1995] (hereinafter referred to as HG) is <D,(E; D") = 
2.8 x 10 s cm 2 s' 1 when averaged over the D + bulge. This flux 
was derived at 500 km. Because the D* densities at low solar 
activity are below the detection range of the OIMS at this 
altitude, we do not have adequate statistics for D* densities to 
determine <b e (E; D*) directly. Instead, we make the reasonable 
assumption that the ratio of the average D + density at low solar 
activity to the one at solar maximum is the same as the ratio 
for H\ in which case e (E; D*) =0.15 for both IT and D + . For 
the D charge exchange escape flux we use the fractionation 
factor/=0.02 [Krasnopolsky, 1985] to obtain d^fCE; D) = 
Rf<t> e (C E; H) = 2.5 x 10 : x 0.02 x 0.9 cm 2 s' 1 =*4.5 x I0 } 
cm * s ', where R is the D/H ratio. Using the scaling relation 
defined above, T(D)rt(D)/[Af(D)MD)] = 0.54, we obtain 
(CE; D) = 2.4 x l O’ cm 2 s 1 and £ (CE; D) = 0.54. Clearly, 
only 0,(E; D*) and <$> e (E; D + ) need be considered in the escape 
of D because the respective escape fluxes due to charge 
exchange are 2 orders of magnitude less. Following a similar 
procedure to that used to derive the source flux for H, a D source 
flux <t> s (D) = 4.3 x 10 5 cm V is obtained in the appendix. 
Using these values, escape efficiencies and globally averaged 
escape fluxes from (9) are shown in Table 4 and compared with 
the flux due to electric field acceleration (HG). Two values of 


Table 4. Escape Efficiencies and Globally Averaged Escape Fluxes GO 5 cm 2 s ! ) 



e 

n 





Process 


(£„ = 0 25) 

(From (8)) 

(£„ = 0.4) 

(From (8)) 

(From HG) 


0.00 

0.75 

1.61 

0.60 

1.29 


E 

0.15 

0.78 

1.68 

0.64 

1.38 

0 56 


1.00 


2.15 


2.15 
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the density ratio have been used for comparison; namely, e n = 
0.25 and t n = 0.4, the average value and the one for Venus year 
2, respectively. 

To produce agreement between the observed HG flux and the 
calculated value, <& 5 (D) would have to be reduced by a factor 
between one third and one half. This could be done by reducing 
the deuterium density at the homopause by about the same 
factor in order to maintain the same density at 220 km with the 
reduced flux. The D/H ratio in the mixed atmosphere would be 
between 0.8 x 10 2 and 1.3 x 10 2 instead of 2.5 x 10' 2 , as 
deduced by Donahue and Hodges [1992]. A ratio as small as 
1.5 x 10' 2 is compatible with one measurement of [HDO]/[H 2 0] 
by near-infrared spectroscopic sounding [deBergh et al., 
1991]. Two other measurements [Donahue and Hodges , 1993; 
Bjoraker et ai , 1992] set the lower limit at 2 x 10 2 . So an 
adjustment of the D/H ratio can at best bring down to about 
1.1 x 10 3 cm' 2 s ' or 0.83 x 10 5 cm' 2 s 1 for z n = 0.25 or 0.4, 
respectively. Furthermore, the value of O^D) used may be 
somewhat larger than the actual one because a fraction of the 
upward flow in the lower thermosphere will be diverted by 
horizontal winds and never flow across the lower boundary of 
the region where the flux balance is calculated. An estimate of 
this effect could be made with a yet to be developed three- 
dimensional circulation model which includes the light 
constituents H and D. The foregoing discussion was prompted 
by the differences in average quantities; however, it should be 


kept in mind that when the scatter in the measurements is 
taken into account, the observed HG flux and the calculated 
flux are then within the variances of one another. 

4.2. Time-Dependent Considerations 

In early 1982 the 10.7-cm solar flux decreased rapidly 
during 4 months from about 220 x 10 22 W m 2 Hz to about 1 20 
x 10' 22 W m 2 Hz (Figure 1). By the time hydrogen and 
deuterium were being observed in the bulge region in mid- 
1992, their densities, as we have seen, were higher than they 
had been during the first three Venus years by a factor of about 
6.5 in the case of hydrogen and a factor of 4 in the case of 
deuterium. Considering only the region of the bulge, the time- 
dependent expression (3) can be used to evaluate the time 
average vertical flux 

«p, + <P</> = 2.5«p s > - A»/A t (11) 

where is the column density of H above the turbopause and 
the ratio A d fA b = 2.5 has been applied. Because the hydrogen 
and deuterium densities were measured for only a small range of 
altitudes near 160 km, it is not possible to evaluate A^/Af 
precisely, nor <p, (r) and <p^(r) separately. Therefore several 
extreme cases will be considered. These demonstrate the 
degree of uncertainty that exists with regard to the 
establishment of the steady state below 160 km. It is possible 
that the densities at and above 160 km are approximately 



H Density in cm * 3 


Figure 7. Altitude profiles of H densities at the subsolar point for various upward fluxes (solid lines). 
Altitude profiles of H densities in the bulge region at solar maximum (short-dashed lipe) and solar minimum 
(long-short-dashed line). 
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Table 5 . Rate of Change of Hydrogen 


Column Density 

(10 7 cm' 1 s' 1 ) 

• 


Case(A/ ) 

AW At 

a(0.8) 

5.6 

b(0.8) 

7.9 

c(0.8) 

10 


’time in years 


those of the steady state and can be used to calculate the 
eventual steady state fluxes, even though the steady state 
downward flux and density below 160 km may not have been 
reached. Thus, the time required for the hydrogen content of 
the bulge to change to its observed value will be taken to be 
either (1) 0.4 year, the time interval during which solar 
activity decreased from its high to low value, or (2) 0.8 year, 
the time from the end of solar maximum to the end of 
observations in 1992. Three possibilities for Afl' will be 
considered. The assumed steady state H bulge profiles are 
shown in Figure 7. These were obtained by extrapolating the 
H densities deduced from measurements at 160 km to 1 10 km, 
where the density is assumed to be the same as in the daytime. 
Above 160 km, diffusive equilibrium is assumed. An extreme 
lower limit to A& is obtained for the first possibility, case a, 
by assuming that the change with time occurs only above 1 60 
km. Case b approximates the effective density in 1992 below 
160 km by extending the H profile above 160 km back to 120 
km, where it meets the solar maximum profile (dotted line), 
and case c assumes that as At -4 0.8 year, a steady state had 
been reached and that AW is the difference between the two 
steady state curves at all altitudes. The results for At = 0.8 year 
are shown in Table 5. The values for a(0.4), b(0.4), and c(0.4) 
are all twice as large as the rates in Table 5. As for < (p 5 >, it can 
be assumed that (p^ increases from d> 5 to the limiting flux 
cm 2 s' 1 as soon as solar activity decreased, 
or that it remains fixed at its solar maximum value. 

The value for «p f + ips obtained by subtracting AW At from 
2.5<<pj> should be larger than the least value that <<p f + <p^> 
can reasonably be expected to have. This can be obtained by 
assuming that <p, decreases from <t> e to 0.15 in 0.4 year and 
then remains at the lower level while <p^ remains fixed at its 
steady state level: 


«Pd> = 2 5 4» s - <t> ( . (12) 

Thus 

«P, + %'min - 2.5*, - 0.425 [1 + l/(p+l)] (13) 

where p , the power of time in the integrand of the time average 
integral, is taken to be 1. Eventually <p rf will become 6.5 
times larger, but as we shall see, such large values in 1992 
cannot be tolerated. In order to maintain consistency in the 
steady state and time-dependent analyses, is taken to be the 
calculated flux in Table 3, not the observed flux in Table 1 . 
These assumptions are consistent with the steady state 
analysis in which d> e was compatible only with the electric 
field escape process for which e =0.15. Thus it is necessary 
that 

2.5<t> Smax - AW/At ^ 2.5d> s - 0.638 (14) 

On the right-hand side of (14), the choice of e 5 makes very 
little difference. The terms in (14) are compared in Table 6. 
This table shows that to meet condition (14) it is necessary 
that steady state conditions for solar minimum were not yet 
established in the lower thermosphere during preentry. Thus 
hydrogen densities below 160 km and <p^ must have remained 
close to their 1991 levels throughout the preentry period to 
allow AW on the left-hand side of (14) and <(p e > on the right 
hand side to be sufficiently small as to satisfy condition (14). 

In the case of deuterium a similar analysis gives the results 
in Table 7, where V is the column density of D above the 
turbopause. When At = 0.4 year, case values are twice as large 
as those in Table 7. <J> 5mar is 6.5 x 10 3 cm' 2 s' 1 . Table 8 
compares terms in the analog of conditions in (14). 

The criterion (14) is easily satisfied for deuterium under a 
wide range of conditions. Thus the results of the time- 
dependent analysis, while not quite as rewarding as those for 
the steady state, are reasonably close to those expected if 
indeed the escape fluxes from Venus, especially for hydrogen, 
are quite small. They can exceed those shown in Table 3 for 
the electric field process E by very little, just as in the case of 
the steady state. 

5. Discussion and Summary 

In this paper, a strong solar cycle variation in the H and D 
densities in the nightside thermosphere of Venus is reported. 
Prominent density enhancements in H of 6.5 times and D of 4 
times are found in the nighttime bulge region as solar activity 
decreases from maximum to near minimum values. Knowing 


Table 6. Comparison of Terms in (14) for Hydrogen (10 7 cm' 2 s ’) 



2 ^Smax- A * lAt 


- 0.638<t> f 


Case 

At = 0.8 yr 

C> 

ii 

o 

4^ 

'C 

Year 

Value 


a 

5.4 

0.0 

1 

3.5 


b 

3.1 

-4.8 

2 

3.3 


c 

1.0 

-9.0 

3 

3.0 





Average 

3.3 
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Table 7. Rate of Change of Deuterium 
Column Density (10* cm 2 s' 1 ) 


Cue(A/«) 

# 

AVI At 

a(0.8) 


2.5 

b(0.8) 


5.0 

c(0.8) 


8.8 


‘time in years 


that the dominant escape fluxes due to the charge separation 
electric field (E) and charge exchange (CE) decrease as solar 
activity declines [Donahue and Hartle, 1992], we propose that 
the observed buildup of H and D is a direct consequence of this 
decrease. 

Making use of the observed H and D densities for both the 
high and low solar activity conditions, the respective flow 
budgets are evaluated with only minimal knowledge about the 
escape processes involved. The H and D escape fluxes so 
derived serve as an independent test of the escape scenario 
described by Donahue and Hartle [1992]. In the following, the 
concept behind the derivation of these escape fluxes is briefly 
summarized and is illustrated in Figure 8. 

Given the observed H and D densities on the dayside, the 
upward fluxes of these constituents, (solar maximum 
conditions) or 0$ (solar minimum conditions) from the 
mesospheric source region can be derived from standard 
diffusion equations. As illustrated in Figure 8, these fluxes 
feed the net horizontal transport (due to winds and exospheric 
return flow) that effectively carries H and D into the nightside, 
where they pile up and attain values that are 2 orders of 
magnitude higher than those on the dayside. Such large 
density enhancements naturally lead to increased downward 
transport having flux magnitudes, <P d or 0^, taken to be 
proportional to the densities with proportionality factors that 
are the same for both high and low solar activity conditions. 
Essentially all escape of H and D is expected to occur from the 
bulge region with flux magnitudes, or 0^, that are 
considered unknown quantities to be derived. The expected 
inverse relations between the H and D bulge densities and 
escape fluxes are indicated in Figure 8. Finally, the H and D 
escape fluxes are uniquely defined when conservation of flow is 
applied simultaneously to both high and low solar activity 
conditions. 


Applying this concept, we developed a simple relationship 
(9) between the planet averaged escape and source fluxes at 
solar maximum ; i.e., , where the escape efficiency 

r| provides the connecting link between parameters at high and 
low solar activity. With knowledge of the source fluxes and 
the solar cycle dependence of the bulge densities, expression 
(9) was used to evaluate the H and D escape fluxes. Source 
fluxes 3>5(H) = 1.45 x 10 7 cm' 2 s'* and 3> 5 (D) = 2.15 x 10 3 
cm 2 s' 1 were derived, where the former is a reevaluation of a 
previous derivation by Paxton et al. [1988]. By spanning the 
range of possible parameters in the escape efficiencies, the 
escape fluxes were found to be in the range 1.23 x 10 7 cm 2 s 1 
£ (H) £ 1.45 x 10 7 cm' 2 s' 1 and 1.61 x 10 3 cm' 2 s' 1 £ <!\,(D) 

£ 2.15 x 10 3 cm' 2 s \ where H and D densities in the respective 
bulges, averaged over the first three Venus years (solar 
maximum), were used. 

It is fortunate that the ranges of escape fluxes obtained by 
analysis of the solar cycle induced bulge changes were found to 
be so narrow that they can be used to analyze the evolution of 
H and D without specifying any particular escape mechanisms. 
Alternatively, these fluxes also provide valuable constraints 
for the study of escape mechanisms. For example, the H escape 
fluxes derived directly from the charge exchange mechanism, 
<E,(CE; H) = 0.9 x 10 7 cm' 2 s’ 1 , and the electric filed 
mechanism, 3> e (E; H) » 1.5 x 10 7 cm 2 s 1 , are just outside the 
bounds derived from analysis of the solar cycle induced bulge 
changes. Thus, considering the variance of the measured 
densities and other possible uncertainties such as the effective 
areas of the bulges where escape takes place, these fluxes and 
their sum <5*(E + CE; H) = 2.4 x i0 7 cm 2 s 1 appear to be 
reasonable estimates. In contrast, the D escape flux due to the 
electric field process, <E,(E; D) = 0.56 x 10 5 cm' 2 s \ is 
outside the experimental range derived from solar cycle 
variations by about a factor of 3 (escape due to CE is not 
important here because the fractionation factor is 0.02). This 
gap can be reduced by a factor of 2 if the experimentally 
possible D/H ratio of 1.5 x 10* 2 instead of 2.5 x 10* 2 is used in 
the derivation of 3>5(D). When the other uncertainties 
mentioned in section 4 are taken into account, the gap should 
be reduced even further. The analysis presented here thus 
provides independent estimates of the H and D escape fluxes 
that are in reasonable agreement with values derived earlier 
[Donahue and Hartle , 1992; Hartle and Grebowsky, 1995]. 

The above analysis assumed that steady state conditions 
applied at the high and low solar activity periods when the 
measurements were made. We have examined the validity of 
this assumption by applying the time- dependent variation 


Table 8, Comparison of Terms in the Analog of (14) for Deuterium (10 6 cm : s' 1 ) 


2.5 ^Smax - AVIAt 2.50> 5 - 0.638<t> e 


Case At = 0.8 yr At - 0.4 yr Year Value 


a 

1.3 

1.04 

1 

0.70 

b 

11 

0.60 

2 

0.89 

c 

0.7 

-0.16 

3 

0.65 




Average 

064 
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Maximum Solar Activity 


Low Solar Activity 


Escape Row (*,) 



Night Day 


Figure 8. Schematic of flow balances in thermosphere 
periods. 


Escape Flow (+ e ) 



described in (3) to the period when solar activity changed from 
its maximum value to the low value at the end of the PVO 
mission. We found that steady state conditions were 
established for H above about 160 km but not yet in the lower 
thermosphere during the preentry period, while they appear to 
be for D throughout the thermosphere. Therefore we conclude 
that the steady state analysis is a valid approximation because 
it dealt with only the observed density variations in the upper 
thermosphere above about 160 km. Consequently, these 
results can be used in an evaluation of the evolution of H and 
D. 


The fractionation factor 


1 <<p e (D)> 
~R <<p e (H)> 


(15) 


is about 0.44, using the results listed in Tables 2 and 4 and R , 
the current D/H ratio of 2.5 x 10' 2 . This large value is 
incompatible with charge exchange for which / is 0.02 
[Krasnopolsky, 1985]. Hence the combined results suggest 
that virtually all of the escape flux determined by this 
procedure is due to the charge separation electric field. In any 
case, taking the H and D fluxes obtained in this way without 
specifying the physical mechanism, we arrive at the solar 
cycle average rates 


«P,(H)> = 3.6 x 10 7 cm 2 s' 1 
<<p tf (D)> = 4 x 10 5 cm' 2 s ' 
/ = 0.44 


when e rt (D) = 0.4. The fluxes averaged over the entire planet 
and the solar cycle are «p/H)> = 7.2 x 10* cm* 2 s * 1 and <<p ( (D)> = 
8 x 10 4 cm' 2 s' 1 . 

The present water vapor mixing ratio is 30 ppm (0.015 m of 
liquid or 8.3 x 10 22 H atoms cm* 2 ). If this is the remnant of 
simple Rayteigh fractionation of an early water reservoir 
where the D/H ratio R Q was terrestrial or 1.6 x 10"* and there are 
no hydrogen sources comparable to «p,(H)> the size of the early 
reservoir was a factor 


r = (R/R c f /(l =8 34jt io 3 (16) 


times as large as the present one. This is enough water to 
cover the planet with 125 m of liquid and is 3.8% of a full 
terrestrial ocean. 

If, on the other hand, hydrogen is in a steady state currently 
with some source (e.g., cometary water and/or outgassed 
interior water) whose D/H ratio is R, = 1.6 x 10* 4 the ratio of 
the primitive reservoir to the present one [Gunvell, 1995; 
Donahue et ai, 1996] is given by 

r = pe fi/l: H -\)j. (17) 

Here p is R/R Q or 157 and 

x H = $/«(>, (H> (18) 

with being the column abundance of hydrogen in the entire 
atmosphere. We note that (17) results when r/x H » 1 in 
Gurwell’s [1995] full equation, which implies that H reaches a 
steady state long before D does. Here, x H is 365 m.y., // x H is 
12.3, and r is 3.8 x 10 4 . The consequence would be 570 m of 
water or 17% of a full terrestrial ocean 4.5 Gyr ago on Venus. 

With such a large fractionation factor, the relaxation time 
for deuterium, x D , is very small and the time required to 
establish a steady state 

T ss = V/ 09) 

is only 830 m.y.. However, the D/H ratio in a source with 
which the present-day hydrogen is in a steady state would have 
to be very large 

R s = R ssf= 11 x 10' 2 . (20) 

If instead the fluxes in Table 1 with the corresponding 
deuterium fluxes are used, we can consider two cases, one i n 
which process CE is added to E, the other in which it is not. In 
the first case,/ = 0.1, and in the second /is 0.15. Then 

x H = (160 - 306) m.y., 
r = 275 or 383, 

calling for an original reservoir 4.1 or 5.7 m deep when there 
is no external source of water. Considering ; = 4.5 Gyr, then tf 
x H would be either 14.7 or 28.1, which calls for a ratio r of 
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Table 9. Three Approaches Used to Determine Escape Fluxes 


Key Parameters 

Solar Cycle Bulge Density Charges 
/ = 0.44 

Electric Field 
/= 0.15 

Electric Field and Charge Exchange 
/ =0.1 

No Source 

Source 

No Source 

Source 

No Source 

Source 

t H (my) 

365 

365 

306 

306 

160 

160 

t/I H 

- 

12.3 

- 

14.7 

• 

28.1 

< v* 

3.6 

3.6 

4.3 

4.3 

8.2 

8.2 

rx 10' 5 

8.3 

38 

0.28 

1.37 

0.38 

2.45 

d (m) 

125 

570 

4.13 

20.6 

5.75 

36,8 


•units of 10 7 cm' 2 s* 1 


1370 to 2450 between the original and present water 
abundance if hydrogen is in a steady state with a source having 
an Earth-like D/H ratio. This amounts to 20.6 or 36.8 m of 
water uniformly distributed on the surface. 

The various scenarios described above for the decay of an 
ancient reservoir of water are conveniently summarized in 
Table 9 in terms of the three approaches used to determine 
escape fluxes; namely, those derived from solar cycle changes 
in the H and D bulges, those obtained from measurements 
related to electric field acceleration, and those produced by 
combining escape fluxes due to charge exchange and electric 
field processes. The table identifies key parameters used in 
describing the decay of an early water reservoir on Venus, 
where time / is 4.5 Gyr, <<p,> is the global average hydrogen 
escape flux, r is the ratio of the initial column of water in the 
atmosphere to the current one, and d is its depth in liquid form 
uniformly distributed over the surface. For each method, the 
size of an ancient water reservoir is determined for two 
possibilities, one with no external source of water and the 
other in which the D/H ratio of the source R s is 1.6 x 10' 4 . In 
all cases, the D/H enrichment p is 157. For each method, the 
amount of water required in an initial reservoir when the 
hydrogen escape flux is in equilibrium with an external source 
is about 5 to 6 times greater than is required for simple decay 
(no source) of the reservoir. This result is expected because the 
primordial water reservoir must be larger in the steady state 
case to overwhelm the low D/H water brought in by an external 
(e.g., cometary) source. Although the hydrogen escape fluxes 
only differ by about a factor of 2, the sizes of the ancient water 
reservoirs differ widely to attain the deuterium enrichment 
observed in Venus’ atmosphere today. This is due to the fact 
that while the fractionation factors / differ by 3 to 4 times, 
these differences are magnified by the exponential dependence 
of an initial water reservoir on /. The differences in the 
fractionation factors are primarily attributed to the differences 
in the respective deuterium escape fluxes. At this time we have 
no way of determining which (p tf (D) is most likely. However, 
future work aimed at minimizing some of the uncertainties 
mentioned above could narrow this gap. For example, since 
the source fluxes and downward fluxes are dependent on global 
circulation, improved estimates of these fluxes and their 
effective areas could be obtained by a three-dimensional, 
global circulation model for the major constituents and H and 


D (e.g., similar to the one developed by Mengel et al [1989] 
for He). This approach would require inclusion of lateral 
transport and escape of H and D, because these processes are 
expected to have important effects on their circulation. Of 
course, experiments to measure the structure and dynamics of H 
and D throughout the thermosphere on a future mission to 
Venus could greatly improve our understanding of escape and 
evolution of the atmosphere. 


Appendix 


Paxton et al. [1988] compared daytime Lyman-a obtained 
in nadir to horizon scans by the Pioneer Venus Orbiter 
Ultraviolet spectrometer (OUVS) during the First three Venus 
years with that predicted from model hydrogen density 
profiles. These profiles were the result of calculating the 
density with a specified upward hydrogen flow <p from the 
expression 


n(z) = n c ( Z ]\ 


j _ | dz!_ 

^ n c (z')[K(z’) + D(z')] 


(Al) 


where K(z) and D(z) are the eddy and molecular diffusion 
coefficients 


/ i T < z) 

n c (z) = n c (z c ) — 
l c 


X exp 


D(z>) . + . K(z>) 


Hj(z') HJz) 


dz 


K(z')+D(z‘) 


(A2) 


Here, H { and H a are the scale heights for hydrogen and the bulk 
atmosphere. Paxton et al. set the reference altitude z c at 200 
km. The emission rates near the horizon are sensitive mainly 
to the density at high altitude, while the emission rates from 
directly below the spacecraft are essentially determined by the 
integrated vertical column density. This, in turn, is set by (p 
and the hydrogen density below the hydrogen homopause. The 
best fit to their data came from a profile with a flux of 7 5 x 10 
cm' 2 s' 1 and a density of 6 x 10 4 cm 3 at 200 km. 

This profile calls for a density of 6 x 10 7 cm 3 at 110 km. 
There is a problem associating a flux as large as 7.5 x 10 7 cm 2 
s' 1 with such a profile because that flux exceeds the limiting 
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flux of hydrogen under the conditions specified. The limiting 
flux <p, is D(100) n,/H a . With VIRA 

D(z) = 8.4 x 10 17 T 0A I n a (z). (A3) 

Here, <p,is 3.3 x 10 7 cm* 2 s*\ when T = 178°K, n a = 8.3 x 10 13 
cm' 3 at 110 km. We have recalculated the hydrogen densities, 
with n = 6 x 10 7 cm* 5 at 1 10 km and K = 1.4 x 10 ,3 /V n a £ 3 x 
10 7 cm 2 s’ 1 . The results for a range of tp from 0 to <p, are shown 
in Figure 7. A density of 6 x 10 4 cm' 3 at 200 km results for 
the profile with <p = 2.3 x 10 7 cm’ 2 s’ 1 . That is a profile very 
similar to that of Paxton et al. with cp = 7.5 x 10 7 cm 2 s’ 1 . In 
this paper we assume that the PVOUVS results for 1978-1981 
are best represented by this profile. Because the average 
hydrogen densities for Venus years 1, 2, and 3 were 7.5 x 10 4 
cm' 3 at 200 km, the densities and fluxes in Figure (7) should be 
multiplied by a factor 1.25 to apply to the analysis in this 
paper. 

The noontime deuterium densities at 220 km during solar 
maximum average 730 cm 3 . If the D/H ratio is 2.5 x 10 2 for 
Venus hydrogen, the deuterium density at 1 10 km will be 1.9 x 
10 6 cm' 3 when the hydrogen density is 7.5 x 10 7 cm 3 . We have 
calculated n a for deuterium from (A2) with the molecular 
diffusion coefficient 

D(z) = 2.3 x 10' 7 I 075 / n a iz) (A4) 

the same as that for H 2 in C0 2 [Hunten, 1973]. The diffusive 
equilibrium density at 220 km according to (Al) is 1.3 x 10 4 
cm' 3 , It follows from (A2) that the observed deuterium density 
of 730 cm 3 at 220 km calls for an upward deuterium flux of 4.3 
x 10 5 cm' 2 s'* near the subsolar point. To reduce that flux 
sufficiently to bring the calculated <P e into agreement with the 
HG flux, n c would have to be reduced by a factor between 0.3 (e 
= 1, year 3) and 0.5 (e = 1.3, year 2). This would call for the 
D/H ratio to be correspondingly reduced. 

Acknowledgment. The authors thank Joel Selekof of Hughes STX 
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Abstract. The hydrogen ion concentrations recently 
observed on Venus, near solar minimum, by the Ion Mass 
Spectrometer on the Pioneer Venus Orbiter in the anti-solar 
sector (22:00-02:00 LST) of the ionosphere are more than an 
order of magnitude less than those previously observed at 
solar maximum. This strong solar cycle variation has a 
profound effect on the escape of hydrogen (and deuterium) 
from Venus; almost all escape occurs during solar maximum. 
After adjustment for solar cycle variation, a planet- averaged 
hydrogen escape flux of 0.6- 1.4 x 10 7 cm" 2 s* 1 is obtained 
along with a large deuterium fractionation factor of 0.1-0.14. 
These results suggest at least two plausible scenarios for the 
evolution of water on Venus: (1) Water vapor on Venus 
may be approaching a steady state if the escape flux is 
balanced by endogenous or exogenous sources of water. The 
source of water must be highly fractionated, with a D/H ratio 
differing by less than an order of magnitude from the present 
ratio of 2.4 x 10" 2, thus precluding low D/H water from 
comets, asteroids or a mantle reservoir. (2) The present day 
D/H ratio of 2.4 x 10* 2 could be established by Rayleigh 
fractionation of an early low D/H water reservoir if the 
escape flux was sufficiently large in earlier times. An early 
water endowment at least 340 times today's abundance, 
equivalent to 4.2 to 14 m of liquid water on the surface, 
would be needed. 

1. Introduction 

The Pioneer Venus Orbiter (PVO) spacecraft, which 
explored the thermosphere of Venus between December, 
1978 and July, 1980, has recently revisited that region. 
During the late summer and fall of 1992 the spacecraft was 
exploring the nighttime atmosphere at altitudes as low as 132 
km. 1978-80 was a time of maximum solar activity. In 
September, 1992, activity was declining toward a minimum. 
Hence the Venus thermosphere and ionosphere was 
examined with the same instruments during times of high and 
low activity. The present paper will contrast the hydrogen 
ion population of the ionosphere between 150 km and 700 
km in the anti-solar sector (22:00 - 02:00 LST) under these 
very different conditions. 

2. Hydrogen Ion Densities at High and Low Solar Activity 

Figure 1 shows examples of hydrogen ion densities 
obtained by the PV Orbiter Ion Mass Spectrometer (OIMS) 
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on orbits 5021 and 5032 (LST 00:47 - 01:57). Ions of mass 2 
are also plotted. These are characteristic of conditions during 
the present epoch. These profiles provide a vivid contrast 
with conditions in 1978-80 when the densities observed 
under similar conditions were about 10 times larger at 200 
km and 30 times larger at 600 km [Bauer et al., 1985]. Even 
though some of these data points will need further scrutiny, 
no fine tuning can change the essential element on which this 
analysis depends. The general level of the densities will not 
change. It is very low. Figure 2 shows average profiles for 
H+ in filled and depleted ionospheres between 22:00 and 
02:00 LST (orbits 4995-5033, 31 cases in all). Averages for 
the two hours before and after midnight are shown separately. 
In Figure 3 the VIRA model H+ densities [Bauer et al., 1985] 
for 150° - 170° SZA are shown along with the current results. 

3. Implications for Hydrogen Escape 

We wish to draw attention to the significance of such large 
variations in H+ for escape of hydrogen and deuterium from 
Venus. Until recently four processes have been recognized 
as contributing currently to escape: classical "Jeans" escape 
(J), charge exchange between low temperature atoms and 
high temperature ions (CE), collisions of hydrogen atoms 
with high speed oxygen atoms produced by dissociative 
recombination of C>2 + (O*), and outflow in the plasma tail 
(T). The first process is virtually inactive in the cold upper 
atmosphere of Venus. Calculations of CE by Hodges and 
Tinsley [1986] and by Rodriguez et al. [1984] give disparate 
results -- 2.4 x 10 7 cm" 2 s’ ^ in the first case and 0.4 - 1 x 10 7 
cm -2 s"l in the second. The latter authors also model a loss 
by O* amounting to 1.2 x 10^ cm -2 s" * that contrasts with a 
recent result of 3.5 x 10^ cm* 2 s'l obtained by Gurwell and 
Yung [1992]. (T) has been discussed by Brace et al. [1987] 
who produced evidence for 0 + escape in 1984 and argued for 
an accompanying — but unobserved — planet wide average 
H + flux of 5 x 10^ ions/cm" 2 s'*. Two remarks need to be 
made concerning the calculations. First, the ionospheric 
models on which they were based were constructed with data 
obtained by the PV orbiter during the solar maximum of 
1978-80 (or the solar minimum of 1984 for T). Second, the 
ionospheric models adapted by the two groups who have 
calculated CE differ strikingly from each other, and neither 
corresponds closely to the VIRA empirical model, which is 
no more than an average of (OIMS and OETP) data obtained 
during three successive passages of penapsis over the night 
hemisphere (Figure 3). Obviously, these calculations need to 
be repeated using H+ densities obtained from PVO 
observations. Allowance should be made for changes with 
solar activity which will reduce the average CE flux by a 
factor of 2 unless compensating changes in the neutral 
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H + Densities in cm * 3 


Fig. 2 Average H + densities in the midnight sector ot 
Venus, 31 cases from Orbits 4995-5032, 9 August - 14 
September, 1992. 

density occur. We note that O will also vary, because radio 
occultation data [Kliore and Mullen, 1988] clearly show that 
there is a large solar cycle modulation in 02 + density. 

While awaiting the results of these exercises, it may be 
possible to reconcile the results already published and to 
adjust them for solar variability by recourse to the PV DIMS 
results. From Figure 3 it is clear that the Tinsley and Hodges 
[1986] ion densities above 250 km are loo high by a tactor ot 
3 f or SZA 150°- 180°. A similar comparison at other zenith 



Fig. 3 H + densities Solar Zenith Angle 150° * 180“. OIMS 
and OETP, 1978-80 (VIRAL OIMS 1992 and models used 
by Hodges and Tinsley (H.T.), [1986] and Rodrigue/ et 
al.(R-P-M). [ 1984]. 


angles shows that, overall, they are effectively high by a 
factor of about 2.7 for the entire nightside. Thus them CE 
probably should be reduced from 2.8 to 1 x 10 7 cm'- s' . 
Likewise the Rodriguez et al. [1984] fluxes are too low by a 
factor of 2.3 and 1.9 for the SZA zones 150°- 1 80° and 120“ - 
150°, in which almost all CE is generated. For their Case 3 
this changes their CE flux from 3.2 x 10 6 cm' 2 s' 1 to 9 x H) 6 
cm' 2 s‘ l during solar maximum. Thus, the results ot the two 
calculations may converge at about 1() 7 cm' 2 s'^ it they are 
performed using the same (observed) H 4 " densities. 

It is interesting to compare these theoretical tluxes with the 
recent determination of the upward flux ot H + in the 
nightside thermosphere of Venus by Hartle and Grebowsky 
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[1992]. Data obtained in the dawn (hydrogen bulge) sector 
between midnight and 2 a.m., 150°- 180° SZA show that H + 
(and D+) ions in this region flow upward between 300 km 
and the ion exobase at 500 km, accelerated by the charge 
separation electric field. This occurs in filled ionospheres 
under quiet solar conditions -- as well as in ionospheric holes. 
The H + flux reaches a local value of 8 x 10 7 cm - s’l and the 
D + flux 3 x 10 5 cm' 2 s' 1 at 500 km, with values larger by a 
factor of 2 — especially for D + -- possible. The midnight 
region was selected because the lateral flow of ions from the 
dayside stagnates there, and the vertical flow is not 
complicated by the effects of day to night transport. There is 
reason to expect that the outward flow detected here will also 
be a component of the vertical flow everywhere on the 
nightside, and be particularly strong in the region of the 
bulge. Allowing for variation in H+ density with SZA in the 
bulge leads to the estimate that this flux would attain a value 
of about 14 x 10 7 cm -2 s' 1 between 150° and 120° and drop 
to about 2.9 x 10 7 cm* 2 between 1 20° and 90°. The average 
bulge flux would be 7.5 x 10 7 cm'- s'K This would 
correspond to a globally averaged flux of 1.5 x 1() 7 cm’ 2 s’ ^ 
if it assumed that almost all of the How occurs in the bulge 
which covers 20 % of the planet's surface. 

These upward flowing ions have two possible destinations. 
One of these is the plasma flow of H+ and D+ accompanying 
the tail ray 0 + flow described by Brace el al. [1987] which 
we have called T. The other is an analog of the telluric polar 
wind in which the upward flowing ions are accelerated to 
escape speed and escape into the solar wind with a flux P. 
Thus a new process (P) would be added to the four already 
identified in contributing to escape. We proceed now to 
estimate the contribution of these five processes to obtain the 
hydrogen escape flux averaged over the planet and over the 
solar cycle. We shall also estimate the fractionation factor 
for each process and for the aggregate of the five. 

4. Escape Fluxes and Fractionation Factors 


The fractionation factor for D escape relative to H is given 


by 


f = 


0D 

a4>n’ 


U) 


where a is the D/H ratio in the mixed atmosphere 
[Krasnopolsky. 1985]. With a = 2.4 x 10* 2 [Donahue et al., 
1982; de Bergh et al., 1991; T.M. Donahue and R.R. Hodges, 
Jr., unpublished 1992] f is 0.17 (and perhaps considerably 
larger) for the thermospheric D+ and H + flows obtained by 
Hartle and Grebowsky [1992]. If the tail ray flow is really a 
plasma flow, the fractionation factor ft is merely the ratio of 

the D+ and H+ densities divided by a. Since the density ratio 
obtained by Hartle and Grebowsky [1992] is 6.2 x 10' 2 , it 
follows that fi = 0.26. If the D+ and H+ ions in the Venus 
wind flow' (P) have equal kinetic energies, as they should if 
they are accelerated by an electric field, fp will be fi/V2 or 
0.15. (We note that the ratio of H + and D+ velocities 
obtained by Hartle and Grebowsky is, in fact, 1.4.) It follow's 
then, that P is 1.2 x 10 7 cm' 2 s’^ and T is 0.3 x 10 7 cm' 2 s'^. 

Combining these flux values with the estimate of 0.9 x 10 7 
cm' 2 s' 1 , arrived at in Section 3 for CE, and 0.35 x 10 7 cm’ 2 


Table 1. Escape Fluxes and Fractionation Factors 


Process 

0 

cm" 2 s'* 

f 

0 

cm' 2 s'l 

f 

f 

CE 

0.9d) 

0 . 02 ( 5 ) 

0.9 

0.02 

0.02 

O* 

0.35(2) 

0.31(2) 

0.35 

o( 4 ) 

0.31 

t 

0.3d) 

0.26 




P 

1.2d) 

0.15(1) 




total 

2.75 

0.14 

1.25 

0.013(5) 

0.1 

average 

1.4 


0.6 




( 1 )This paper; ( 2 )Gurwell and Yung, [1992]; 
(•^Krasnopolsky, [1985]; ^McElroy et al., [1982]; 
(SlHunten et al., [1982].. 


s'l for O* calculated by Gurwell and Yung [1992] gives a 
total flux of 2.75 x 10 7 cm* 2 s' 1 at solar maximum, or a solar 
cycled averaged flux 0 of 1.4 x 10 7 cm' 2 s‘^. Taking the 
values derived for fp and ft together with 0.02 for fee and 
0.31 for f G * [Gurwell and Yung, 1992] yields a total 

fractionation factor of f =0.14. Here we have always taken 
a to be 2.4 x 10‘ 2 or 150 times terrestrial. 

In Table 1 we display the contributions of various 
mechanisms to the escape flux and fractionation factors that 
have been associated with these processes. The first two 
columns give the values suggested in this paper. The next 
three columns show the effect of neglecting tail ray and 
Venus wind losses but allowing for the large increase in f 0 * 
proposed by Gurwell and Yung [1992], From this table it is 
apparent that only charge exchange, and Jeans escape now 
seem to discriminate severely against deuterium escape. No 
matter w'hat combination of processes is considered the 
fractionation factor is 0.1 or larger, and the effective 
hvdrogen escape flux is severely reduced by solar cycle 
effects. We shall now discuss the implications of these 
changes for the evolution of water on Venus. 

5. Evolution of Venus Water 

Recently there has been a convergence of some 
measurements of the water vapor mixing ratio in the lower 
atmosphere toward 30 ppm [Pollack et al., 1992; Donahue 
and Hodges, unpublished results 1992]. However, some 
older credible analyses that call for much higher mixing 
ratios — up to 200 ppm — remain unexplained. Thus there is 
a possibility that the water vapor concentration is variable in 
space and time and that possibility should still be allowed for 
in a discussion of the water vapor budget of Venus — present 
and past. The rest of the hydrogen inventory is virtually 
unbounded except by aeronomic constraints. Here we shall 
consider a possible range of 30 to 100 ppm H 2 O equivalent. 
The time required to exhaust such reservoirs by the escape 
flux discussed here is given by 

T = H/0, (2) 

where H is 8.3 x 10 22 atoms enr 2 if the water vapor mixing 
ratio is 30 ppm. For a flux of 1.4 x 10 2 cm' 2 s'* t lies 
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between 190 and 630 My. The characteristic time for decay 
of the D/H signature as D/H approaches a steady state value 
in the presence of sources and sinks of water, is 

t» = t/? (3) 

[Krasnopolsky, 1985]. This is between 1.4 and 4.5 Gy. The 
order of magnitude increase in f resulting from the 
reevaluation of fo* and the large values of ft and fp, 
therefore, means that it is plausible that the present day water 
vapor is well on its way toward a steady state if sufficiently 
strong sources -- endogenous or exogenous — exist. 
However, the same high efficiency for deuterium escape 
which, on the one hand, makes x$$ short, on the other hand, 
makes the ultimate D/H ratio, given by 

(D/H)« = ou/f , (4) 

[Krasnopolsky, 1985] only 7 times the source value, a s . 
Because water in comets very likely has a D/H ratio in the 
neighborhood of a few times 10“^, this would appear to 
preclude comets [Grinspoon and Lewis, 1988] as important 
sources and to require mechanisms for producing a highly 
fractionated source of water from a mantle reservoir, given 
the present very high D/H ratio of 2.4 x 10'^. 

On the other hand it would be possible for the presently 
large value of D/H to be established by simple Rayleigh 
fractionation of an early water reservoir with a terrestrial like 
D/H ratio if ij) should increase adequately with H. If ct/a 0 is 
150 and f is 0.14 an ancient water supply 340 times larger 
than the present atmospheric inventory would be required. 
This is equivalent to 0.14 to 0.47 percent of a full terrestrial 
ocean. To exhaust the associated hydrogen, 0 would have to 
increase appropriately with increasing water abundance 
[Kumar et al., 1985]. 

6. Conclusions 

Solar activity has a profound effect on hydrogen escape 
from Venus; the fractionation factor for deuterium and 
hydrogen escape is very large and our understanding of the 
history of water on Venus is very confused. Virtually the 
only scenario that appears to be excluded by our present 
understanding of hydrogen and deuterium escape processes is 
one that involves a steady state in which H and D outflow are 
balanced by an input of low D/H water from comets, 
asteroids or a mantle reservoir. However, it would be 
difficult to underestimate the quality of our present 
understanding either of the processes or the present hydrogen 
inventory. Developments are occurring, it remains to be seen 
whether or not they represent progress. 

Conclusions, similar to those reached here regarding the 
significance of the large fractionation factor generated by 
process O* and T, were reached, prior to the preparation of 
this paper, by D. H. Grinspoon [Unpublished manuscript; 
International Colloquium on Venus, Aug. 10-12, Pasadena, 
CA.] 
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